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FOREV'ORD 

This  final  technical  report  entitled,  "A  Determinis¬ 
tic  Approach  to  the  Prediction  of  Tcleseismic  Ground  Motion 
From  Nuclear  Explosions,"  is  submitted  by  Systems,  Science 
and  Software  (S3)  to  the  Advanced  Research  Projects  Agency 
(AREA)  and  to  the  Defense  Nuclear  Agency  (DMA).  The  report 
presents  the  results  from  a  twelve-month  effort  to  determine 
the  sensitivity  of  M  and  m^  to  variations  in  explosion 
yield,  depth  of  burial,  near  source  material  properties  and 
earth  structure.  The  work  was  performed  under  Contract 
Number  DNA  001-73-C-0167 .  Col.  J.  T.  Jones  was  the  ARPA 
Program  Manager  and  Mr.  Clifton  B.  McFarland  was  the  DNA 
Project  Scientist. 

Dr.  J.  Theodore  Cherry  was  the  S3  Project  Manager  for 
the  study.  Professors  Charles  B.  Archambeau  and  David  G. 
Ilarkr ider  of  the  California  Institute  of  Technology  were 
consultants  on  the  project.  Their  interaction  with  Dr. 
Thomas  C.  Bache,  Dr.  Gerald  A.  Frazier  and  Mr.  Kenneth  G. 
Hamilton  have  made  the  teleseismic  codes  viable  prediction 
tools. 
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I.  INTRODUCTION 

This  report  presents  theoretical  predictions  or  tele- 
seismic.  ground  motion  from  underground  nuclear  explosions. 

The  theoretical  predictions,  which  are  presented  in  the  form 
of  synthetic  seismograms,  are  obtained  from  a  deterministic 
computer  model  which  merges  nonlinear  shock  code  calculations 
in  the  near-field  with  seismic  code  calculations  of  far-field 
stress  wave  propagation  through  a  realistic  earth  model.  The 
application  of  this  computer  model  to  study  dependence  of 
teleseismic  ground  motion  on  explosion  yield,  depth  of  burial, 
near-source  material  properties,  tectonic  stress  release  and 
the  earth  structure  between  source  and  receiver  is  the  sub¬ 
ject  of  this  report. 

Assuming  spherical  symmetry,  the  explosion  representa¬ 
tion  in  terms  of  an  equivalent  elastic  source  (the  reduced 
displacement  potential,  RDP)  is  calculated  using  a  one¬ 
dimensional,  Lagrangian  finite-difference  code  (SKIPPER). 

This  shock  wave  code  requires,  as  input,  a  description  of 
the  near  source  material  response  to  arbitrary  stress  load¬ 
ing.  The  shock  wave  from  the  high  pressure  cavity  which  con¬ 
tains  the  rock  vaporized  by  the  explosion  is  then  propagated 
to  a  radius  at  which  the  response  is  linear  elastic.  The  RDP 
is  then  obtained  from  the  calculated  displacement  or  velo¬ 
city  history  of  a  particle  in  the  elastic  region. 

The  teleseismic  bedy  and  surface  wave  codes  accept 
the  SKIPPER  calculated  RDP  and  propagate  it  through  a  given 
heterogeneous  earth  model.  The  body  wave  code  uses  ray 
theory  in  a  radially  stratified  earth  model  while  the  sur¬ 
face  wave  code  assumes  a  flat  earth  model  and  sums  the 
vertically  trapped  modes  of  propagation.  The  shock  wave 
code  and  its  interaction  with  the  body  and  surface  wave 
codes  has  been  described  by  Cherry,  et  al.^  The  details 
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of  the  surface  wave  calculations  are  given  by  Harkrider . ^ 2 ^ 

Tlie  ray  theory  used  in  the  body  wave  calculation  nay  be 
found  in  Pullen. ^ 

In  Section  II  of  this  report  the  calculation  of  the 
equivalent  elastic  source  (RDP)  for  explosions  in  tuff, 
granite  and  sandstone  are  described.  A  study  of  the  sur¬ 
face  waves  emanating  from  explosions  in  these  three  media 
is  presented  in  Section  III.  In  Section  IV  the  correspond¬ 
ing  body  wave  calculations  are  presented. 

A  major  goal  of  this  investigation  has  been  to  deter¬ 
mine  the  dependence  of  body  and  surface  wave  magnitudes,  m^ 
and  Ms ,  on  the  important  parameters  characterizing  the  ex¬ 
plosions  which  are  listed  in  the  first  paragraph.  Since 
theoretical  seismograms  at  teleseismic  distances  are  now 
routinely  obtained,  and  are  picked  directly  from  the 
seismograms  in  a  standard  manner  and  compared  to  observations. 

In  this  way  the  model  parameters  controlling  mb  and  M.  are 
identified. 

For  a  given  earth  structure,  is  found  to  be  directly 
proportional  to  log^  <K°°)  ,  the  logarithm  of  the  steady  state 
value  of  the  RDP.  On  the  other  hand,  body  wave  magnitude  is 
found  to  be  strongly  influenced  by  depth  of  burial  and  near 
source  compressional  velocity,  as  well  as  i^O).  The  shape  of 
the  source  spectrum  has  no  influence  on  M  and  only  a  minor 
effect  on  m^ .  Also,  while  tectonic  stress  release  may  add 
considerable  complexity  to  the  P  wave  train,  it  has  little 
effect  on  m^ .  For  both  body  and  surface  waves,  the  importance 
of  an  accurate  earth  structure  model  is  emphasized. 

The  deterministic  approach  used  in  this  research,  fea¬ 
turing  a  synthesis  of  results  from  rock  mechanics,  linear  and 
nonlinear  continuum  mechanics,  computational  physics  and  seis¬ 
mology,  now  permits  a  detailed  analysis  of  the  factors  controlling 
the  teleseismic  signature  of  a  nuclear  explosion. 
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II.  CALCULATED  SOURCE  FUNCTIONS 


2.1  THE  REDUCED  DISPLACEMENT  POTENTIAL 

SKIPPER^’  ^  is  a  Lagrangian,  finite-difference 
code  that  numerically  simulates  a  propagating  stress  (shock) 
wave  in  one  space  dimension.  The  code  is  capable  of  carrying 
the  stress  field  into  the  small  displacement;  elastic  regime 
and  yet  flexible  enough  to  permit  appropriate  material  response 
formulations  in  the  large  displacement  nonlinear  regime. 

For  a  spherically  symmetric  explosion  source  the 
reduced  displacement  potential,  iJj,  is  given  by 


U(t,r) 


where 


<L_  ^(t)  =  }P_  + 

3r  r  ^2  rot  » 


(2.1) 


U(t,r)  =  radial  particle  displacement  at  t  and 
r  in  the  elastic  region 

a  =  the  compressional  velocity  in  the  earth 
medium 

t  =  time 


r  =  radial  position 
i  =  t  -  r/a  =  retarded  time 

ip  =  dip/ dr  =  the  reduced  velocity  potential  (RVP) 

Since  ip  is  only  a  function  of  t  in  the  elastic 
region,  its  wave  form  does  not  change  with  changing  r. 
Therefore,  Eq.  (2.1)  may  be  used  to  determine  ip  at  any 
radius  in  the  elastic  tone.  A  subroutine  exists  within  the 
SKIPPER  code  that  numerically  integrates  Eq.  (2.1)  at  various 
specified  radii  from  the  center  of  the  explosion.  The  fact 
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that  ip  should  be  invariant  with  1  he  radius  provides  an 
important  verification  that  the  specified  radii  are  in  the 
elastic  region. 

The  following  equations  show  the  relation  between  the 
RDP  (and  its  derivatives)  and  particle  velocity  (u) ,  radial 
stress  (Trr)  ,  tangential  stress  (TQe)  and  the  hydrodynamic 
component  of  stress  (P): 


u  = 


T 

rr 


-P 


a 


2P 


nr 


+  46 


(Jl  +  Jl_\ 

'r3  r2a  ' 


(2.2) 

(2.3) 


(2.4) 

(2.5) 

where 

k  =  the  bulk  modulus, 

U  =  the  shear  modulus, 
p  =  density, 


00 


k  -  2/3y 


■)  P  +  2p(^-  +  1-  ) 

'  '  r 3  r2a/ 


p  =  - 


a 

6 


=  Jr'+ 

*  P  ’ 

~  =  shear  wave  velocity. 


These  equations  show  that  only  the  hydrodynamic  component 
of  stress  will  attenuate  as  r"1  over  the  entire  elastic 
region.  Also,  from  Eq.  (2.1)  the  steady-state  value  of 
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the  RDP  is  given  by 


1H00)  =  r  2U  (°°,  r)  (26) 

where  U(»,r)  is  the  steady-state  (static)  displacement  at 
the  radius,  r. 

Both  the  body  wave  and  surface  wave  teleseismic  codes 
accept  the  equivalent  source  function  in  the  frequency  do¬ 
main.  Since  the  RDP,  i p,  has  a  finite  steady-state  value,  it 
Fourier  transform  is  difficult  to  evaluate  numerically.  The 
RVP,  has  a  steady-state  value  of  $,(»)  =  0,  so  that  its 
Fourier  transform  is  bounded  at  zero  frequency.  If  $  is 
the  Fourier  transform  of  the  RVP,  then  at  zero  frequency 


(v  =  0)  = 


^(T)dT  =  f  3T  dT  • 

*0 


Therefore,  since  ip (0 )  =  0, 


■HO)  =  *(")  (2.7) 

All  our  calculated  source  functions  will  be  presented  in  the 

^  format.  Equation  (2.7)  relates  an  important  time  domain 

parameter,  iK“) ,  to  the  value  of  the  RVP  transform  at  zero 
frequency . 


We  will  find  it  convenient  to  scale  the  source  func¬ 
tion,  calculated  from  SKIPPER,  to  different  explosion  yields 
Table  2.1  gives  the  appropriate  cube  root  scaling  for  each  o: 
the  listed  variables.  This  scaling  follows  from  the  assump¬ 
tion  that  cavity  volume  (L3)  is  directly  proportional  to 
yield  (W) .  Therefore  L  ^  W1/3.  Since  the  wave  velocity 
in  the  medium  is  independent  of  yield,  then  T  v 


he  have  found  that  the  scaling  given  in  Table  2.1 
is  in  fact  applicable  to  SKIPPER  calculations  having  the  same 
properties  but  differing  greatly  in  the  explosion  yield. 
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TUFF 


The  technique  used  in  SKIPPER  to  obtain  the  pressure 
component  (P)  of  the  equation  of  state  has  been  presented 
in  detail  in  previous  reports^  ’  and  only  a  brief  summary 
will  be  given  here. 

The  code  initially  assumes  a  pressure  equilibrium  mix- 
ing  of  grain  density  tuff  and  water,  where  the  mix  is  obtained 
by  specifying  f,  the  mass  fraction  of  water.  Figures  2.1a  and 
2.1b  show  the  Ilugoniot  and  a  feur  release  isentropes  for  these 
two  constituents.  Figure  2.2  shows  the  result  of  mixing  the 
constituents  of  Fig.  2.1  using  f  =  0.17. 

Partial  saturation  is  developed  by  specifying  both 

the  initial  void  fraction,  $  ,  and  the  rate  of  irreversible 

o 

void  collapse  during  pressure  loading.  Figure  2.3  shows  the 

loading  and  release  P-v  curves  that  result  from  adding  five 

percent  air-filled  porosity  (<J>  =  0.05)  to  the  mix  of  Fig.  2.2. 

0 

Also  shown  in  the  figure  are  the  experimentally  determined 
unloading  data  for  Unit  3  tuff  at  the  NTS  Diamond  Dust  site.^ 
The  pressure  state  labeled  Pc  in  the  figure  represents 
the  state  at  which  total  irreversible  pore  collapse  occurs, 
while  Pg  gi\es  the  (elastic)  pressure  at  which  irrever¬ 
sible  pore  collapse  begins. 

We  have  assumed  that  material  strength  depends  on  the 
pressure  (P)  and  internal  energy  (e)  and  is  governed  by  a  re¬ 
lation  of  the  form 

SijSji  ^  T  Y2CP.e),  (2.8) 


where  T^  -  <5  „  P  is  the  stress  deviator  tensor  and 


V(f,e)  =pP  2-t  1  -  f- 

m  m 


P  <  P 


P  >  P  , 
m 


(2.9) 


e  >  e 

m 


i 


t  o 


re  (kbar) 


SKIPPER  Loading 


-Loading  and  release  P-v  curves  for  par 
*  =0.05). 
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For  a  given  stress  state  in  which  two  of  the  principal 
stresses  are  equa  1 »  corresponds  to  the  maximum  allowable 

stress  difference  that  the  material  is  able  to  support.  Re¬ 
ported  values  of  Ym  for  the  tuff  from  the  Diamond  Mine 
and  Diamond  Dust  sites  range  between  0.4  and  2.0  kbar.^ 
he  have  used  the  range  as  a  guide  for  the  input  variable  Y  . 

For  spherical  symmetry,  Eq.  (2.8)  becomes 


where 
tion.  When 


S  is  the  deviatoric  stress 
rr  2 


5rr  |  exceeds  j  Y,  at 

culation  then  we  set 
2 


S 


rr 


3  V  (sign  S  J 


component  in  the  radial  direc- 
some  point  in  the  SKIPPLR  cal- 

(2.11) 


This  equation  of  state  formulation  has  been  used  in 
order  to  determine  the  important  material  properties  that 
control  the  source  spectrum.  ^  It  has  also  been  used  to  pro¬ 
vide  an  excellent  pre-shot  prediction  of  the  ground  motion  for 
the  Mine  Dust  HE  test.  ^ 


Most  of  the  nuclear  explosions  in  NTS  tuff,  having  yields 

large  enough  to  give  an  adequate  sampling  of  ground  motion  at 

teleseismic  distances,  were  below  the  water  table.  Figure  2.4 

1/3 

shows  DOB  versus  W  '  for  those  events  having  unclassified 
[781 

yields.  *  J  BILBY  was  located  in  Area  3  at  NTS.  Pahute 
Mesa  was  the  location  of  the  remainder  of  the  shots  shown  in 
the  figure. 

Since  the  degree  of  saturation  depends  on  permeability, 
there  is  no  assurance  that  the  rock  environment  below  the  water 
table  is  completely  saturated.  Hence, it  is  instructive  to  show 
the  effect  of  air-filled  porosity  of  the  RVP.  Appropriate 
source  calculations  have  been  performed  by  Cherry,  et  al . , ^ 
and  are  listed  as  T-l,  T-2,  T-3,  T-4  in  Table  2.2.  Calculations 


16 


. . rnmmm 


R-  2  034 


T-5  and  T-6  are  new  and  are  intended  to  show  the  effect  of  over¬ 
burden  pressure  on  the  RVP.  All  the  results  listed  in  the  table 
are  for  a  hypothetical  source  yield  of  0.02  kt. 

TABLE  2 . 2 


SUMMARY  OF  SKIPPER  CALCULATIONS  IN  TUFF 


In  Table  2.2 


<{>o  =  volume  fraction  of  air-filled  porosity 

Ym  =  maximum  material  strength 
=  initial  density 
Pq  =  overburden  pressure 
^cav  =  calculated  final  cavity  radius 

^el  =  calculated  radius  at  which  the  stress  wave  enters 
the  elastic  region,  i.e.,  the  radial  extent  of 
the  "shatter"  zone. 

A 

’KO)  =  zero  frequency  value  of  RVP;  also  equals  late 
time  value  of  RDP ,  i^(oo) 

DOB  =  depth  of  burial  =  P  (g  p)  "  1 

Ip  foo)  0  j 

ml"mi  =  log10  ^  (oo)  =  approximate  (assumes  a  flat  spectrum 
and  no  ^surface  reflection)  change  in  Mg  and  mb 
if  material  is  changed  from  T-l  to  T-i. 

Parameters  common  to  all  calculations  were 
f  =  mass  fraction  of  water  =  0.17 
a  =  compression  velocity  =  2.4  km/sec 
Pg  =  elastic  pressure  =  0.075  kbar 
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P  =  crush  pressure  ^  0.5  kbar 

c 

U  =  shear  modulus  =  40  kbar 

p  =  pressure  at  which  material  sustains  maximum 

m 

strength  (Y^)  =  0.6  kbar. 

Figure  2.5  compares  the  calculated  source  functions 
|,*,|  for  T-l,  T-3  and  T-4.  While  the  calculations  were  mn 
using  a  source  yield  of  0.02  kt ,  the  scaled  frequency  axis 
corresponding  to  a  yield  of  1000  kt  is  also  given  in  the 
figure  (see  Table  1.1).  Increasing  the  porosity  from  zero 
■ t-i)  t0  o.05  (T-3)  decreases  the  RVP  by  a  factor  of  3.6 
for  frequencies  less  than  1  Hz  and  yields  up  to  1000  kt.  If 
in  addition,  the  maximum  material  strength,  Ym,  is  increased 
by  a  factor  of  2.0  (T-l  to  T-4),  the  RVP  is  reduced  by  a 
factor  of  5.5. 

/x 

The  effect  of  the  various  source  functions  |  if  |  of  Fig 
2.5  on  body  wave  and  surface  wave  magnitude  are  given  in 
Table  2.2  in  the  m^-nw  column.  From  the  table  we  find  that 
increasing  the  dry  porosity  from  zero  (T-l)  to  0.02  (T-2) 
causes  a  magnitude  decrease  of  0.31.  The  effect  of  air; 
filled  porosity  on  if(0)  is  shown  in  Fig.  2.6.  Since  if(0) 
scales  as  the  device  yield  (Table  1.1),  then  1  kt  in  T-l  is 
equivalent  to  3.7  kt  in  T-3. 

Finally  the  effect  of  increased  depth  of  burial  on 
if(0)  is  available  from  calculations  T-4,  T-5  and  T-6.  In¬ 
creasing  the  DOB  from  2110  ft  (P  =  116  bars)  to  4500  ft 
(P  =  247  bars)  causes  a  magnitude  decrease  of  only  0.07 
units  (0.74  -  0.81).  As  shown  in  Section  IV,  the  major 
DOB  effect  on  m^  is  caused  by  the  interference  between 
P  and  pP. 
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Frequency  (HZ)  (0 . 0 2  kT) 


1.0  2.0  3.0 

Frequency  (HZ)  (lOOO  kT ) 


Comparison  of  source  functions  for  0.02  kt  from 
Calculations  (T-l),  (T-3)  and  (T-4). 
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2 . 3  SANDSTONE 

Rulison  and  Rio  Blanco  were  nuclear  explosive,  gas 
stimulation  experiments  in  sandstone.  The  Rulison  experiment 
was  located  in  the  Rulison  basin  of  Colorado,  where  a  single 
40  kt  device  was  detonated  at  a  depth  of  8442  ft.  Rio  Blanco 
was  in  the  Piceance  Basin  where  three  30  kt  devices  were  simul¬ 
taneously  detonated  in  a  vertical  emplacement  hole  at  depths 
of  5840,  6230  and  6690  feet. 

An  equation  of  state  for  sandstone  was  developed  from 

r  q  i 

quasistatic  test  data  published  by  Schock,  et  al ■ ,  1  J  and 
high  pressure  sandstone  ilugoniot  data  published  by  Terhune.^^ 
Figure  2.7  shows  pressure  versus  specific  volume  developed  from 
the  quasistatic  tests.  Figure  2.8  merges  the  static  data  in 
the  P-v  plane  writh  the  high  pressure  Ilugoniot  data.  The  maxi¬ 
mum  stress  difference  allowable  under  shock  loading  was  taken 

to  be  Y  =  4.3  kbar.  This  was  assumed  to  occur  at  P  10 
m  m 

kbar.  The  material  strength  relation  is  given  by  Eq .  (2.9). 

TABLE  2.3 

SUMMARY  OF  SKIPPER  CALCULATIONS  IN  SANDSTONE 


Calculation 
Symbo 1 

a 

km/sec 

Y  (kbar) 
m 

Pm(kbar) 

P  (gm/cc) 
0 

P  (kbar) 

0 

Rcav<m> 

Rel(») 

DOB  (ft) 

C  (0) 

(m!) 

S-l 

3.7 

4.3 

10 

2.52 

0.5 

2.3 

16 

6300 

S-2 

3.7 

4.3 

10 

2.52 

0.6 

2.2 

14 

8400 

1  .s 

The  results  of  Table  2.3  are  for  a  hypothetical  source 
yield  of  0.02  kt  and  may  be  compared  directly  with  those  listed 
for  tuff  in  Table  2.2.  S-l  corresponds  to  a  DOB  appropriate 
for  Rio  Blanco  while  S-2  was  run  at  the  Rulison  depth. 

Figures  2.9  and  2.10  show  the  calculated  source  func¬ 
tions  for  S-l  and  S-2.  These  sandstone  source  functions  are 
similar  to  the  T-3  RVP  for  tuff  shown  in  Fig.  2.5.  Given  the 
same  source  function  for  different  materials,  body  wave  ampli¬ 
tudes  will  depend  only  on  the  compress ional  velocity,  a,  in 
Eq.  (2.1). 
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0.34  0.35  0.36  0.37  0.38  0.39  0.40 

v  (cc/gm) 


Fig.  2 . 7- -Pressure  (P)  versus  specific  volume  (v)  for  sand- rg, 
stone  from  the  static  test  data  of  Schock,  £ t  a 1 . 
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°-26  0.28  0.30  0.32  0.34  0.36  0.38 

v  (cc/gn) 


Pressure  (P)  versus  specific  volume  (v)  for  sand 
stone,  merging  the  low  pressure  data  of  Fig.  2.7 
with  the  high  pressure  Hugoniot  data  reported  bv 
Terhune . [10] 
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The  results  given  in  Tables  2.2  and  2.3  a] so  show  that 

the  final  cavity  radius  (P.^)  is  not  simply  related  to  $(0) 

or,  equivalently,  iK«0  .  For  example,  S-l  and  T-5  have  almost 

identical  values  of  tH°°)  but  quite  different  final  cavity 

radii.  It  is  therefore  important  that  Eq  .  (2.6)  only  be 

applied  to  the  elastic  region  defined  bv  r  >  R 

—  '  el  ’ 

2.4  GRANITE 

Granite  source  functions  have  been  calculated  based 
on  the  NTS  granite  compressibility  data  of  Stephens  and 
Lilley  and  the  strength  data  of  Heard. [12J  Figure  2.11 
is  taken  from  Heard  and  shows  the  effect  of  water  content 
(pore  pressure)  on  the  strength  of  Hardhat  granodiorite . 

The  dry  granite  is  three  times  stronger  than  wet  granite 

in  which  the  pore  (fluid)  pressure  equals  the  confining 
pressure . 

Two  calculations  were  run  at  hypothetical  yields  of 
0.02  kt  changing  only  the  material  strength.  The  calculation 
G-l  used  the  dry  strength  data  of  Fig.  2.11  and  G-2  used  the 
wet  strength  data  shown  in  the  figure.  The  source  func¬ 
tions  for  these  two  calculations  are  shown  in  Figs.  2.12 

and  2.13.  The  results  of  the  calculations  are  summarized 
in  Table  2.4. 


TABLE  2.4 

summary  of  skipper  calculations  for  granite 


cnlculat ion 
Symbol 

a 

km/sec 

Yn(kbar) 

Pm(kbnr) 

0,  (gm/cc) 

T  (kbar) 

0 

RcavW 

RclC"0 

POB  (ft) 

G  - 1 

5.6 

20 

22 

2.67 

0.120 

1,82 

33.6 

1500 

G-  2 

5 . 6 

6.67 

22 

2.67 

0.120 

2.51 

65.9 

1500 

The  overburden  pressure,  P^,  was  chosen  to  correspond 
to  the  DOB  of  Piledriver,  a  60  kT  event  in  granite  at  NTS. 
Free  field"  ground  motion  from  this  event  was  processed  by 
i.illiam  R.  Perret  at  Sandia  Laboratory.  The  reduced 
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Fig.  2 . 11 -- Failure  envelopes  for  Ilardhat  granodiorite  tested 
in  compression,  25°C,  at  a  strain  rate  of  10'4/sec,  from 
Heard  [12] . 
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displacement  potentials  that  he  obtained,  scaled  to  0.02  kt 
at  Station  2403  (r  =  2000  ft),  Station  1503  (r  =  2803  ft) 
and  Station  2803  (r  =  2803  ft),  are  shown  in  Figs.  2.14,  2.15 
and  2.16.  The  calculated  RDP  from  G-l  and  G-2  are  given  in 
Figs.  2.17  end  2.18. 

The  agreement  between  the  calculated  and  observed  re¬ 
duced  displacement  potentials  is  poor.  Even  though  >jj(«) 
from  G-2  compares  reasonably  well  with  iii(°°)  from  Stations 
-403  and  2803  (8.2  m3  versus  10.8  m3  and  9  m3),  the  time  at 
which  is  attained  is  at  least  an  order  of  magnitude 

different  between  the  calculated  source  function  G-2  and 
that  observed  at  the  two  stations. 

A  possible  explanation  for  the  discrepancy  is  that  the 
observed  ground  motion  at  the  three  stations  is  not  really 
free  field,  i.e.,  additional  arrivals  have  obscured  the  true 
nature  of  the  equivalent  source  from  the  explosion.  For 
example,  the  free  surface  reflection  should  arrive  at  Station 
2403  at  approximately  0.015  sec,  where  this  time  has  been 
scaled  to  0.02  kt.  Referring  to  Fig.  2.14,  this  arrival 
occurs  during  the  excursion  of  the  RDP  to  its  peak  value  of 
16.1  m  .  Release  of  tectonic  stress  would  be  another  possi¬ 
bility  and  may  explain  the  azimuthal  variations  of  the  RDP 
from  Station  1503  and  2803. 

If  the  observed  ground  motion  at  the  three  stations 
is  indeed  contaminated  by  additional  events,  then  none  of 
the  source  functions  from  these  stations  will  explain  the 
teleseismic  observations  from  Piledriver  and  we  will  have 
to  depend  on  calculated  source  functions  to  explain  the  tele- 
seismic  data. 
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2 . 14- -Reduced  displacement  potential  (RDP)  at  Station 
for  Piledriver,  scaled  to  0.02  kt . 
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III.  SURFACE  WAVES  FROM  UNDERGROUND  EXPLOSIONS 

3.1  INTRODUCTION 

The  surface  waves  radiated  by  underground  explosions 
in  general  include  both  Rayleigh  and  Love  waves.  A  computa¬ 
tional  method  for  generating  synthetic  seismograms  for  these 
surface  waves  has  been  developed  by  Harkrider1  J  and  made 
operational  at  S3.  The  explosion  source  functions  described 
in  Section  II  may  be  input  to  the  surface  wave  computer 
codes  and  the  resulting  teleseismic  surface  waves  are  then 
generated. 

Treating  the  explosion  as  a  spherically  symmetric 
source  (e.g.,  the  SKIPPER  calculation),  only  surface  waves 
of  the  Rayleigh  type  are  generated.  The  surface  wave  sig¬ 
nature  is  then  dependent  on  the  steady  state  value  of  the 
reduced  displacement  potential  O(“0)  and  the  velocity  pro¬ 
file  of  the  crust  and  upper  mantle  between  source  and  re¬ 
ceiver.  The  Rayleigh  wave  seismic  record  for  this  case  is 
only  weakly  dependent  on  other  parameters  such  as  the  burial 
depth  of  the  explosion. 

Many  studies  of  surface  waves  from  underground  ex¬ 
plosions  [e.g.,  13,  14]  have  concluded  that  some  mechanism 
other  than  a  spherically  symmetric  source  was  operational. 
Indications  of  this  secondary  phenomenon  have  included  strong 
Love  wave  radiation  and  an  azimuthal  dependence  of  Rayleigh 
wave  amplitudes.  A  superposition  of  a  double -couple  (quad- 
rupole)  source  on  the  point  source  (monopole)  due  to  the 
explosion  has  been  successful  in  matching  these  observa¬ 
tions.^  a  theory  of  tectonic  stress  release  from  ex¬ 

plosions  in  prestressed  media  has  been  formulated  by 
r  1 7 1 

Archambeau1  1  to  calculate  the  double  -  couple  contribution 
to  the  source.  This  theory  has  been  included  in  our  sur¬ 
face  wave  computational  capability. 
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leleseismic  surface  waves  and  the  measurement  of  sur¬ 
face  wave  magnitude,  fl^ ,  are  sensitive  to  both  the  explo¬ 
sion  and  the  tectonic  stress  release  contributions  to  the 
source.  Many  parameters  enter  into  a  specification  of 
tectonic  release  and  a  detailed  study  of  this  contribution 
v*as  not  conducted  for  the  surface  waves.  Our  primary  attention 
is  directed  to  isolating  the  pure  explosion  component  and  deter¬ 
mining  its  effect  on  teleseismic  surface  waves. 

3 . 2  TELESE I  SMI C  SURFACE  WAVE  CODES 

The  technique  for  propagating  surface  waves  in  a 
layered  earth  medium  is  described  by  ilarkrider.  ^  The 
earth  structure  is  modeled  using  homogeneous  isotropic 
layers  and  the  response  of  the  multilayered  system  is  cal¬ 
culated  using  linear  wave  propagation  theory  for  the  verti¬ 
cally  trapped  modes  of  propagation.  The  eigenfunction  can 
be  calculated  for  any  depth  in  the  model. 

The  source  is  specified  by  the  spherical  wave  multi¬ 
pole  coefficients  describing  an  equivalent  elastic  source. 

Foi  explosions  with  tectonic  release,  only  the  monopole  and 
quadrupole  terms  are  required.  The  spherical  wave  coeffi¬ 
cients  are  then  converted  to  the  corresponding  cylindrical 
wave  expansion  coefficients  associated  with  the  set  of 
potentials  employed  in  the  surface  wave  programs.  These 
coefficients  then  serve  as  excitation  factors  in  the  Love 
and  Rayleigh  wave  programs  and  combine  with  the  eigen¬ 
functions  for  the  layered  system  to  give  predicted  surface 
wave  spectra  at  teleseismic  distances. 

The  synthesis  code  uses  the  eigenvalues  and  eigen¬ 
functions  appropriate  to  the  explosion  depth  to  calculate 
the  spectra  at  selected  distances  from  the  source.  Anelas- 
tic  absorption  along  the  propagation  path  is  computed  using 
values  of  observed  wave  decay  constants.  The  effect  of  the 
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earth's  curvature  on  geometric  spreading  is  also  included. 

Finally,  synthetic  seismograms  are  generated  by  a 
numerical  transformation  of  the  far  field  spectrum  into  the 
time  domain.  Suitable  seismograph  amplitude  and  phase 
response  data  may  be  included  to  construct  a  synthetic 
seismogram  for  a  given  instrument. 

3-3-  RAYLEIGH  WAVES  FROM  UNDERGROUND  EXPLOSIONS 

The  source  functions  calculated  by  the  SKIPPER  code 
(Section  II)  may  be  input  to  the  surface  wave  codes  to  ob¬ 
tain  teleseismic  Rayleigh  waves.  The  important  parameters 
controlling  this  synthesis  are  the  limiting  values  of  the 
RDP  the  earth  structure  and,  weakly,  the  depth  of 

burial.  The  effect  of  these  parameters  is  studied  in  this 
section . 

3.3.1  Earth  Structure 

The  influence  of  changing  the  velocity  profile  of 
the  crust  and  upper  mantle  is  indicated  by  performing  several 
calculations  using  different  models.  Velocity-depth  profiles 
for  the  three  earth  structure  models  studied  are  given  in 
Fig.  3.1.  All  three  models  are  intended  for  the  structure 
in  the  vicinity  of  NTS  and  are  therefore  rather  similar. 

Synthetic  seismograms  at  4000  km  from  the  three  earth 
structure  models  are  presented  in  Fig.  3.2.  These  seismo¬ 
grams  include  the  effect  of  an  LRSM  Long  Period  seismometer. 
The  relative  magnification  of  this  instrument  is  shown  in 
Fig.  3.3.  The  seismograms  are  appropriate  to  the  Piledriver 
event  which  was  a  61  kt  explosion  at  a  burial  depth  of  463 
meters  in  NTS  granite.  The  SKIPPER  source  function  G-l  was 
used  (see  Section  2.4). 
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(a)  Earth  structure  CIT  109. 


to.  3. 2- -Rayleigh  wave  synthetic  seismograms  for  the  Piledrivei 
event  using  three  earth  structure  models.  SKIPPER  source  G-l 
epicen  ral  distance  4000  km,  depth  of  burial  463  meters,  ERsV 
Long  Period  Seismograph.  The  20  second  phase  on  the  radial 
motion  record  is  indicated  (t). 
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The  earth  structure  affects  both  the  arrival  time  and 


amplitude  of  the  various  frequency  components.  Surface  vaves 
from  different  events  are  often  compared  through  a  single 


number,  the  surface  wave  magnitude,  M  . 

[T81 

Gutenberg-Richter  formula  for  is L  1 


The  commonly  used 


M  =  log  A,,  +  1.656  log  A  +  1.818  , 
s  H 


(3.1) 


where  A^  is  the  maximum  zero-to-peak  amplitude  in  microns 
of  the  horizontal  component  of  the  fundamental  mode  Rayleigh 
wave  having  a  period  of  about  20  seconds,  and  A  is  the 
epicentral  distance  in  degrees.  This  formula  is  to  be  ap¬ 
plied  only  for  15  <  A  <  130. 

Several  other  definitions  of  Mg  are  commonly  used. 
In  this  report  we  will  refer  to  a  modification  of  (3.1)  to 
deal  with  maximum  amplitudes  different  from  20  seconds 
which  is, 

Ms  =  log  |  — j  +  1.656  logA  +  3.119  .  ( 


Far  T  =  20,  (3.2)  is  identical  to  (3.1). 

The  LRSM  instrument  for  which  the  seismograms  of 
Fig.  3.3  were  prepared  is  normalized  to  unity  at  25  seconds. 
The  instrument  effect  can  be  deleted  by  applying  an  instru¬ 
ment  correction  to  the  calculated  values,  using  the 

nominal  magnification  of  Fig.  3.3.  The  instrument  corrected 

M  values,  which  represent  true  ground  motion,  will  be 
S  Q 

denoted  M  . 

In  Table  3.1  the  M  values  for  the  records  of  Fig. 
3.2  are  summarized.  From  the  table  we  see  that  the  earth 
structure  can  have  a  noticeable  effect  on  Mg,  even  for  these 
three  similar  structure  models.  A  more  significant  variation 
in  M  results  from  choice  of  the  phase  of  maximum  ampli¬ 
tude,  %he  Airy  phase,  for  Mg  measurements,  rather  than  the 


44 


R- 2034 


TABLE  3.1 

Ms  VALUES  Fnn.  TIT  SEISMOGRAMS  HF  FIGURE  3.3 


|  Structure 

Period 

Ah(m) 

M 

s 

mg 

S 

CIT  109 

20  sec 

14  sec 

0.125 

0.329 

3.49 

4.07 

3.52 

4.34 

NTS  Granite/H.M. 

20  sec 

14  sec 

0.105 

0.235 

3.42 

3.92 

3.45 

4.19 

Fallon  to  Ruth/H.M. 

20  sec 

14  sec 

0.038 

0.210 

3.34 

3.87 

3.37 

4.14 

amplitude  at  20  seconds.  Values  of  M  from  the  two  are 
really  not  comparable. 


A  sampling  of  observed  Mg  values  for  Piledriver  from 
LRSM  data[19]  is  given  in  Table  3.2.  These  values  were  cal¬ 
culated  by  using  the  LRSM  observed  A/T  values  for  the  LQ 
phase  recorded  on  LPR  instruments  in  Eq.  (3.2),  The  observed 
Ms  values  are  somewhat  greater  than  the  computed  values,  which  is 
expected  since  only  explosion  effects  were  considered  in  the 
calculations  leading  to  Table  3.1.  A  strong  tectonic  release 
component  was  observed  for  the  Piledriver  event  ancl  addi¬ 
tion  of  this  contribution  would  raise  our  M  values. 

Further,  the  source  function  G-2  has  a  which  is  3,9 

times  greater  than  the  iH00)  for  G-l,  As  will  be  verified 
in  subsequent  sections,  use  of  G-2  will  increase  all  M  , 

Ms  values  in  Table  3.3  by  0.59. 


At  present  we  are  limited  to  laterally  homogeneous 
structure  models.  Propagation  of  surface  waves  through 
continental -oceanic  boundaries  or  into  regions  of  quite  dif¬ 
ferent  crustal  structure  can,  of  course,  result  in  average 
velocity-depth  profiles  quite  different  from  those  of  Fig,  3.1 
and  have  a  significant  effect  on  M  measurements. 
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TABLE  3.2 

OBSERVED  M  VALUES  FOR  PILEDRIVER 


5 • 3 • 2  Depth  of  Burial 

The  Rayleigh  wave  calculation  leading  to  the  seismo- 
g  of  Fig.  3.2a  is  repeated  at  a  burial  depth  of  1200  m. 
The  two  seismograms  demonstrating  depth  of  burial  effects' 
are  shown  in  Fig.  3.4.  Comparing  Fig.  3.4  to  Fig.  3.2a, 
the  shape  of  the  two  is  seen  to  be  identical,  with  the  only 
ciange  being  a  small  decrease  in  amplitude  with  increasing 
depth.  The  corresponding  Mg  values,  which  are  virtually 
identical,  are  compared  in  Table  3.3. 


TABLE  3.3 

COMPARISON  OF  M  FOR  TWO  DEPTHS  OF  BURIAL 


Rr2034 


3  •  4  RAYLEIGH  WAVES  FROM  SELECTED  EVENTS 

The  teleseismic  surface  waves  from  several  under¬ 
ground  explosions  m  tuff,  sandstone  and  granite  media  will 
now  be  studied.  Since  surface  waves  are  dependent  on  the 
low  frequency  behavior  of  the  source,  and  since  the  cal¬ 
culated  source  functions  are  flat  in  the  low  frequency 

(f  <  0.2  Hz)  portion  of  the  spectrum,  the  surface  wave  ampli¬ 
tudes  are  expected  to  scale  with  the  steady  state  value  of 
the  RDP,  ijj(eo). 

The  basic  parameters  characterizing  several  represen¬ 
tative  underground  explosion  events  in  the  three  source  media 
are  listed  in  Table  3.4.  The  value  of  is  calculated 

for  the  indicated  yield  by  scaling  the  SKIPPER  value  which 
was  calculated  for  a  nominal  yield  of  0.02  kt.  For  example, 
from  Table  2.2  the  steady  state  RDP  for  T-4  is  1.35  m3 . 


TABLE  3.4 

DECf!RIPTIOI!  OF  SELECTED  EXPLOSION  EVENTS 


Event 

Source 

Medium 

Yield 

Depth 

SKIPPER 

Source 

♦  O) 

(m3) 

Bilby 

Tuff 

^200  kt 

700  m 

T-4 

1 .35x10“ 

Rio 

Blanco 

Sand- 
s  tone 

3x30  kT 

1830  m 

S-l 

9 . 0x10  3 

Rul ison 

Sand¬ 

stone 

40  kT 

2573  m 

S  -  2 

3 . 6x10  3 

Pile- 

driver 

Granite 

61  kT 

463  m 

G-l 

6.77xl03 

Knicker¬ 

bocker 

Tuff 

75  kT 

630  m 

T-4 

5 . 06 2x 1 0 3 

Greeley 

Tuff 

825  kT 

1214  m 

T-4 

5.57x10“ 

Benham 

Tuff 

1100  kT 

1402  m 

T-4 

7.43x10“ 
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Using  the  scaling  rule  of  Table  2.1,  the  for  Bilhy  is 

=  Ply  =  13,500  m3  . 

200  U<  Z  0.02 


A  sampling  of  the  vertical  and  horizontal  Rayleigh 
wave  synthetic  seismograms  for  Bilby  and  Rio  Blanco  are 
given  in  Figs.  3.5  and  3.6.  Similar  seismograms  were  con¬ 
structed  for  Rulison.  The  earth  structure  CIT  109  was  used 
for  these  calculations .  The  comparable  seismograms,  using 
the  same  earth  structure,  for  the  Piledriver  event  were 
previously  shown  in  Fig.  3.1a.  Theoretical  magnitudes  for 
these  four  events  are  summarized  in  Table  3.5. 

The  peak-to-peak  amplitude,  2A„,  from  Table  3.5  is 
plotted  vs  iH00)  in  Fig.  3.7.  The  line  of  unit  slope 
through  each  data  set  is  also  shown  on  the  figure.  Recall 
from  Section  3.3.2  that  the  effect  of  increasing  depth  of 
burial  is  a  small  decrease  in  Rayleigh  wave  amplitude. 
Correlating  the  data  of  Fig.  3.7  with  the  depths  for  the 
events  from  Table  3.4,  we  see  that  the  deviation  from  the 
line  of  unit  slope  is  consistent  with  expected  depth  of 
burial  effects. 


Since  Rayleigh  wave  amplitudes  scale 
'P  O)  ,  keeping  depth  of  burial  fixed,  we  can 
estimates  for  the  tuff,  events  Knickerbocker, 
Benham  by  scaling  the  data  of  Tables  3.3  and 
suits  are  summarized  in  Table  3.6. 


directly  with 
obtain  M 

s 

Greeley  and 
3.5.  The  re- 
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-Synthetic  Rayleigh  wave  seismograms  for 
Source  T-4,  DOB  700  meters,  earth  struc 


Fig.  3.5  (Cont . 
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e  seismograms  for  the  Rio  Blanco  event.  Viel 
ve rage  HOB  1S30  meters,  earth  structure  model 
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Fig.  3.6  (Cont. 
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*  Airy  Phase  at  3000  kn 
A  Airy  Phase  at  4000  km 

■  20  sec  Phase  at  3000  km 

•  20  sec  Phase  at  4000  lcn 


Rulison  |  |  Bilby 

Rio  Blanco 
Piledriver 

iK°°) 


ig.  3.7--Peak  to  peak  amplitude  of  the  radial  Rayleigh  wave, 
Aji,  for  selected  events  versus  the  steady  state  value  of  RPP, 
(°°)  •  Lines  of  unit  slope  are  shown  with  each  data  set. 


tabu:  3.5 


PREDICTED  MAGNITUDES  FOR  SELECTED  EVENTS 


Event 

Distance 

(km) 

Period 

(sec) 

2Au(y) 

M 

s 

MC 

S 

Bilby 

3000 

20 

0.72 

3.75 

3.78 

3000 

14 

1.7 

4.27 

4  .  54 

4000 

20 

0.49 

3.8 

3.82 

4000 

14 

1.24 

4.34 

4.61 

Rio  Blanco 

3000 

20 

0.50 

3.59 

3.62 

3000 

14 

1.10 

4.08 

4.35 

4000 

20 

0.32 

3.60 

3.63 

4000 

14 

0.78 

4.14 

4.41 

Rulison 

3000 

20 

0.18 

3.14 

3.17 

3000 

14 

0.40 

3.64 

3.91 

4000 

20 

0.137 

3.23 

3.26 

4000 

14 

0.29 

3.71 

3.98 

Piledriver 

4000 

20 

0.25 

3.49 

3.52 

4000 

- 4 

14 

— 

0,66 

4.07 

4.34 

TABLE  3.6 


PREDICTED  SURFACE  WAVE  RESULTS  FOR  KNICKERBOCKER,  GREELEY 

AND  BENHAM  AT  4000  KM 


Period 

G 

Event 

(sec) 

Ajj  (u) 

Ms 

Mu 

s 

Knickerbocker 

20 

0.093 

3,36 

3.39 

14 

0.233 

3.92 

4.19 

Greeley 

20 

0.987 

4.39 

4.42 

14 

2.54 

4 . 96 

5.23 

Benham 

20 

1.32 

4.52 

4.55 

l 

14 

3.39 

5.08 

5.35 
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COMPARISON  OF  OBSERVED  AMP  COMPUTED 


Treating  each  event  as  a  spherically  symmetric  ex¬ 
plosion  via  the  SKIPPER  calculation,  theoretical  values  for 
Ms  have  been  computed  for  several  representative  events. 

A  good  comparison  of  these  results  with  observed  values 
requires  further  effort  to  collect  observed  data  along  with 
information  on  exactly  how  those  data  were  reduced.  However, 
some  Ms  data  are  on  hand  and.  we  shall  present  the  compari¬ 
son,  recognizing  the  possibility  that  inconsistent  data  re¬ 
duction  exists  over  the  set. 

A  bar  graph  comparison  of  observed  and  computed  M 
is  presented  in  Fig.  3.8.  The  observed  value  for  Renham  is 
from  the  1/WSSN  reports,  that  for  Rulison,  Knickerbocker, 
Greeley  and  Bilby  from  Harkrider  (private  communication) 
and  that  for  Piledriver  from  Table  3.2.  Since  different 
methods  may  have  been  used  in  computing  M  ,  the  observed 
values  are  probably  not  consistent.  The  lower  (left)  points 
on  the  computed  value  bars  correspond  to  M  based  on  the 
20-second  phase. 

Examining  Fig.  3.8,  we  see  that  our  Mg  estimates  are 
consistently  low  compared  to  the  observations.  This  is  not 
surprising  since  the  effect  of  the  anomalous  shear  wave  com¬ 
ponent  (tectonic  release),  present  to  some  degree  in  every 
event,  is  to  raise  Mg .  The  addition  of  a  tectonic  release 
quadrupo le  source  to  the  explosion  source  will  then  increase 
our  computed  Mg .  The  computational  capability  to  accomplish 
this  task  is  available  to  us,  but  has  not  yet  been  applied 
for  surface  waves.  Also,  from  Section  II,  note  that  the  source 
function  1-4,  which  has  been  used  for  our  tuff  calculations,  ha 
the  smallest  iK00)  of  the  tuff  source  functions  presented. 
Selection  of  one  of  the  greater  source  functions,  such  as 
T-3,  would  scale  the  Mg  for  the  tuff  events  accordingly. 
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-Comparison  of  observed  and  computed  M  values  for  selected  events. 


For  example,  source  T-3  results  in  a  new  M  for  the  tuff 
events  which  is  0.2  units  higlier  than  those  shown  in  Fig. 
3.8.  A  similar  upward  scaling  of  M  by  changing  source 
functions  for  the  Piledriver  calculation  was  mentioned  in 
Section  3.3.1. 

In  summary,  selection  of  a  larger  explosion  source 
function,  inclusion  of  tectonic  release  and  modifying  the 
assumed  earth  structure  are  all  effective  ways  to  increase 
the  predicted  values  of  M  .  Before  any  of  these  are  in¬ 
voked,  it  is  important  that  an  M  data  set  be  obtained 
that  is  self-consistent. 
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IV.  TELESEISMIC  BODY  WAVES  FROM  UNDERGROUND  NUCLEAR  EXPLOSIONS 

4.1  INTRODUCTION 

The  teleseismic  signature  of  body  waves  excited  by  an 
underground  nuclear  explosion  is  affected  by  the  various  para¬ 
meters  which  characterize  the  source,  the  receiver,  and  the 
eartli  structure  through  which  the  explosion  generated  stress 
wave  is  propagated.  A  computational  model  which  embodies  the 
main  features  of  these  parameters  is  now  available. 

Synthetic  seismograms  which  provide  accurate  predictions 
of  m.  and  the  basic  shape  of  the  first  few  seconds  of  the 
ground  motion  record  at  selected  teleseismic  distances  are 
now  routinely  generated  by  computer  codes  operational  at  S3. 

The  application  of  these  codes  to  study  the  dependence  of  m. 

(and  the  time  domain  character  of  body  wave  arrivals)  on  the 
important  parameters  affecting  teleseismic  ground  motion  is 
the  subject  of  this  section. 

4.2  THE  BODY  WAVE  CODES 

The  program  for  computing  teleseismic  body  waves  is 
composed  of  three  subprograms  called  SOURCE,  TFRFNC  and  TIMSER. 
The  role  these  subprograms  play  in  the  body  wave  program  is 
the  subject  of  the  following  discussion. 

Description  of  the  Source 

The  explosion  source  (which  may  include  tectonic  stress 

relaxation  effects)  is  described  by  the  displacement  spectra 

(P,  SI!  and  SV)  at  a  specified  radial  distance  r  (in  the 

elastic  regime),  over  a  range  of  azimuthal  angles  0  and  over 

a  range  of  take-off  angles  <j>.  Equivalently,  the  source  may 

be  defined  by  the  coefficients  of  the  spherical  wave  function 

expansion  of  the  outgoing  waves.  For  a  spherically  symmetric 

explosion  plus  tectonic  release,  only  the  monopole  and  quad- 

1171 

rupole  compoents  are  present  in  the  expansion.  1 
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The  reduced  displacement  potential  (RDP)  transform, 
'K^),  is  obtained  from  an  appropriate  SKIPPER  calculation 
as  discussed  in  Section  II.  The  P  wave  displacement  spect¬ 
rum  at  a  radius  r  due  to  a  spherically  symmetric  explosion 
may  be  shown  to  be  related  to  the  RDP  transform  by,^ 


(4.1) 


where  the  superscript  Pg  denotes  the  P  wave  from  the  ex¬ 
plosion,  the  subscripts  r,  0,  <j>  are  spherical  coordinates, 
and  kp  e  w/vp  is  the  wave  number  for  congressional  waves. 
For  a  spherically  symmetric  explosion,  only  P  waves  are  re¬ 
leased.  The  program  SOURCE  reads  J(a>)  and  calculates  the 
displacement  spectrum  uj>  (r,w)  due  t0  the  explosion. 

The  theory  of  the  dynamic  relaxation  of  tectonic 
stress  due  to  the  creation  of  a  cavity  and  surrounding  zone 
of  fractured  material  in  a  prestressed  geologic  formation  is 
given  by  Archambeau.  In  this  theory  the  tectonic  com¬ 

ponent  of  the  source  is  specified  in  terms  of  a  quadrupole 
point  source.  Such  an  equivalent  point  source  representation 
has  been  shown  by  Archambeau  and  others  (e.g.,  Refs.  [13-17]) 
to  explain  the  anomalous  Rayleigh  wave  radiation  pattern  and 
Love  wave  radiation  from  nuclear  explosions.  Calculation  of 
the  displacement  spectra  at  a  radius  r  from^such  a  source 
is  performed  by  SOURCE  yielding  the  spectra  U?t,  USt ,  where 
1  =  r,  6,  <f>  and  the  superscripts  Pt  and  St  indicate  the 
P  and  S  waves  due  to  the  tectonic  release. 


Near-Source  Free  Surface  Reflection 

Having  specified  the  equivalent  elastic  source,  in  a 
suitable  format,  the  remainder  of  the  calculation  is  carried 
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out  in  the  linear  regime.  Therefore,  individual  frequency 
components  may  be  separately  propagated  through  the  earth 
model  and  superposed  at  the  receiver. 

The  program  SOURCE  also  computes  the  reflected  and 
converted  (S  to  P  and  P  to  S)  waves  at  the  free  surface  for 
take-off  angles  of  interest.  Therefore,  if  U?e  and  A  (r,w) 
are  the  P  wave  displacement  spectra  due  to  the  explo¬ 
sion  and  tectonic  release  which  emanate  from  the  source  re¬ 
gion  at  a  take-off  angle  0p  and  aximuth  <j>p,  the  program 
SOURCE  will  compute  the  additional  displacements  UP?e , 

U  t,  from  the  free  surface  reflections  which  also  leave 
the  source  region  along  this  ray  (Fig.  4.1).  An  analogous 
calculation  is  carried  out  for  S  waves. 

For  computation  of  the  near-source  free  surface  re¬ 
flection,  SOURCE  considers  the  source- to-surface  material  to 
be  homogeneous  and  characterized  by  a  separately  input  atten¬ 
uation  parameter  Q  and  wave  speeds  vss  and  vSs.  The 
latter  velocities  are  commonly  selected^  be  tlAsource 
region  values  from  the  SKIPPER  calculation. 

Earth  Structure 

As  just  described,  SOURCE  calculates  the  displacement 
spectrum  for  each  of  the  various  phases  emanating  from  the 
source  region  at  selected  aximuths  and  take-off  angles.  In 
propagating  these  displacements  to  teleseismic  distances,  the 
effect  of  the  earth's  structure  enters  via  multiplication  by 
the  complex  transfer  function,  Ap  (a>) ,  which  includes  the 
radial  spreading,  attenuation  and  dispersion  due  to  the  path 
from  source  to  receiver.  While  several  methods  for  obtain¬ 
ing  such  transfer  functions  are  available,  the  results  of 
this  report  were  obtained  using  the  program  TFRFNC. 
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This  program  uses  ray  theory^“0^  to  calculate  t lie 
transfer  functions  for  ray  paths  through  a  radially  strati¬ 
fied  earth  model  over  a  range  of  take-off  angles.  In 
addition  to  geometric  spreading  effects,  the  wave  amplitudes 
are  attenuated  as 

U  =  Ue'wT/(2  Qp),  (4.2 

where  Q  is  t lie  dissipation  constant  for  the  ray  path 
corresponding  to  tlie  take-off  angle  0  .  A  necessary  conse¬ 
quence  of  tlie  medium  absorption  (4.2)  is  dispersion.  The 
relation  between  the  two  is  discussed  bv  Futterman and 

r  2  ?  1 

Strick.  Using  tlie  formulation  of  Strick,  dispersion  is 

included  by  multiplying  U  of  (4.2)  by 


exp 


TTQp  £n 


(4.3) 


where  u>^  is  selected  as  a  high  frequency  cutoff  parameter 
and  T  is  the  travel  time  associated  with  waves  of  frequency 
Thu s ,  if  Ap  is  the  multiplication  factor  due  to  radial 
spreading,  the  total  transfer  function  for  tlie  ray  of  take¬ 
off  angle  0p  is 


•ojT 


Ap(w)  =  Ap  exp  i  -  jp  £n  (wh/w) 


(4.4) 


Ground  Motion  at  the  Receiver 


Multiplying  by  the  transfer  function  A  (w)  in  SOURCh 
the  spectra  of  the  various  direct,  reflected  and  converted 
displacement  waves  due  to  the  explosion  (and,  if  required, 
tectonic  release)  at  points  of  interest  in  the  teleseismic 
field  are  obtained.  By  applying  suitable  correction  factors 
(see  Appendix  2  of  Ref.  [1]),  these  may  be  converted  into 
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the-  spectra  of  vertical,  radial  (from  ground  cero)  and  tangen¬ 
tial  ground  motion  at  the  receiver.  True  ground  motion  Tor 

each  of  the  phases  results  from  a  lourier  synthesis  into  the 
time  domain. 

The  program  TIMS HR  applies  a  fast  Fourier  transform 
routine  to  the  spectra  of  the  individual  phases  arriving  at 
the  receiver.  If  the  seismometer  transfer  function  is  in¬ 
cluded  prior  to  the  Fourier  transformation,  the  ground  motion 
record  appropriate  to  a  given  instrument  is  obtained. 

Ihe  seismogram  for  each  mantle  arrival  is  obtained  by 
superposing  the  time  records  for  the  phases  associated  with 
a  particular  take-off  angle  (and  ray  path).  Since  several 
take-off  angles  (and  mantle  arrivals)  may  be  associated  with 
a  given  point  in  the  teleseismic  field,  rather  complex  syn¬ 
thetic  seismograms  are  obtained. 

4.3  INVESTIGATION  OF  THE  PARAMETERS  AFFECTING  m 

— - — — _ _ _ b 

4-3-l  Illustrative  Calculations 

Before  proceeding  with  a  detailed  examination  of  the 
various  parameters  affecting  teleseismic  ground  motion  re¬ 
cords,  let  us  describe  a  few  typical  calculations. 

Source  Effects 

Figure  4.2  shows  a  typical  equivalent  source  repre¬ 
sentation,  this  one  appropriate  to  the  Handley  event  which 
was  1000  kt  detonated  at  a  depth  of  1.2  km  in  tuff.  The 
SKIPPER  calculation  for  this  event  is  denoted  T-4  and  was 
scaled  to  the  Handley  yield  using  cube  root  scaling  (Table 
2.1).  The  total  displacement  pulse  and  spectrum  are  shown 
as  well  as  the  same  quantities  including  only  far  field  terms 
(terms  decreasing  as  r" 1  in  Eq .  (4.11)).  Only  far  field 
spectra  and  wave  forms  are  important  for  teleseismic  body 
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unvc  phases.  Figure  4.5  shows  the  theoretical  far  field 
spectrum  and  elastic  wave  form  predicted  on  the  basis  of  a 
different  tuff  material  model  which  is  denoted  T-l.  These 
two  source  material  models,  T-4  and  T-l,  arc  representative 
of  the  variations  of  material  properties  for  tuff. 

The  radiated  wave  form  and  signal  spectrum  due  to  tec¬ 
tonic  release  depend  on  the  effective  size  of  the  prestressed 
region  around  the  induced  zone  of  failure.  This  region  is 
characterized  by  an  effective  radius  R  This  radius  can 
he  taken  to  he  very  large,  effectively  infinite,  or  finite 
and  on  the  order  of,  say,  ten  times  the  shatter  zone  dimen¬ 
sion.  These  two  choices  result  in  different  radiated  seismic 
waves  as  illustrated  in  Fig.  4.4,  which  also  shows  how  this 
difference  is  reflected  in  the  amplitude  .spectrum  of  the 
pulse.  From  the  figure  we  sec  that  a  finite  R_,  reflecting 
a  concentration  of  the  initial  prestress  field  or  strong 
spatial  variation,  results  in  a  peaked  far  field  spectrum 

and  a  distinctive  pulse  shape.  On  the  other  hand,  for  R 

9  s 

infinite  the  far  field  spectrum  is  flat  and  the  signal  pulse 
non -oscillatory.  It  is  difficult  to  assess  which  of  these 
conditions  is  most  common  or  likely  and  we  leave  the  issue 
with  the  comment  that  the  choice  of  R  lias  a  nearly  un¬ 
detectable  effect  on  the  synthetic  seismograms  for  short 
period  P  waves  (the  instrument  filters  out  most  of  the 
difference  and  the  rest  is  obscured  by  simultaneously  arriv¬ 
ing  phases ) . 

Tclcscismic  Spectra 

Multiplying  the  source  spectra  (e.g.,  Fig.  4.5a)  by 
the  transfer  function  characterizing  the  earth,  the  source 
field  is  propagated  to  a  tclcscismic  distance.  In  Fig.  4.5 
the  spectra  due  to  the  Handley  explosion  is  shown  at  a 
distance  of  4066  km.  The  earth  model  used  is  denoted  CIT 
109  Low  Q  (sec  Section  4.5.2).  Both  the  direct  P  wave  and 
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c  "rcleasc\vi>t  h°’’  VT  P)  ™v0  f°™s  and  spectra  fro,, 
P  Thn  ,  1  "s  ^ rclaxat ion  zone  radius)  finite  or 

;c  ,;vcs  represent  true  ground  displacement  with 
tLc  offset  removed.  The  source  dimension  is  approx i- 
with  a  prestress  (shear)  of  65  bars.  The  scales 
ropriate  for  a  distance  of  about  4,000  km  with  the 
tincture  and  0  accounted  for. 
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(a)  Direct  P  wave  from  the  explosion. 


PrMutnc*  (i:*) 


(b)  Reflected  P  wave  fr on  the  explosion 


««  (It.) 

(Cl  ^lc,f,oirCCt  and  reflected  P  kave5  due  to  th«  (d>  Sun  Of  J  i  roc  I ,  rrflrctrd  and  converted  , 
,  #  to  all  sources  (explosion  f»  tectonic  re  lea 

fig.  4 .  5- -Tlicoret  ical  spectra  of  explosion  generate  i  conprcs 
sional  arrivals  from  the  Handley  event  at  a  distance  of  'non 
kn  and  azimuth  of  30°.  CIT  109  Low  0  model,  explosion  model 
T-l,  source  depth  1.2  kn. 
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pP  reflected  from  the  free  surface  are  included.  Figure 
1.5c  shows  the  spectrum  of  P  and  pP  combined.  The  spectral 
peaking  and  scalloped  nature  of  the  high  frequency  portion 
of  the  composite  spectrum  result  from  interference  between 
the  diioct  and  reflected  w a v c s .  In  Fig.  4.5d  the  total 
spectrum  due  to  both  the  explosion  and  tectonic  components 
of  the  source  are  shown.  This  spectrum  is  somewhat  more 

peaked  and  complex  than  Fig.  4.5c  due  to  the  tectonic  release 
effects . 

Figure  4.6  shows  the  individual  spectra  for  the  body 
wave  phases  resulting  from  tectonic  release  (with  R  finite), 
along  with  the  complete  spectrum  of  the  composite  signal  from 
both  explosive  and  tectonic  source  components  (same  as  Fig. 

4 . 5 d ) .  In  this  case  the  largest  effect  is  due  to  the  con¬ 

verted  S  to  P  wave,  since  tectonic  release  is  primarily  a 
source  of  shear  wave  energy.  The  sPT  phase  can  be  as  large 
or  larger  than  the  direct  P  wave  from  the  explosion  itself, 
depending  on  the  stress  field  orientation  and  magnitude  as 
well  as  on  the  size  of  the  shatter  zone  created  by  the  shock 
wave . 

Seismograph  Modulation 

The  spectra  of  Figs.  4.5  and  4.6  relate  to  actual  ground 
displacement  spectra  at  a  point  on  the  free  surface  of  the 
earth.  In  order  to  obtain  a  meaningful  comparison  with  ob¬ 
servations,  this  signal  must  be  modulated  by  the  transfer 
function  of  the  detector.  Since  we  are  interested  in  com¬ 
paring  predictions  of  the  telescismic  field  to  LRSM  station 
recordings,  the  transfer  functions  for  the  short  period 
seismometers  employed  at  these  stations  will  be  used  in  this 
report  Figure  4.7  shows  the  amplitude  and  phase  response 
of  the  LRSM  SP  instrument.  The  gain  has  been  normalized  to 
unity  at  1  Hz.  Multiplying  the  spectra  of  Figs.  4.5  and  4.6 
by  this  transfer  function  and  performing  an  inverse  Fourier 
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(a)  Direct  5'  wave  from  tectonic  release. 


(c)  Converted  S  to  P  wave  from  teutonic  release. 


(b)  Reflected  P  wave  fron  tectonic  release. 


(d)  Sun  of  direct,  reflected  and  cor  rtc.l  ;  ‘ 

to  all  sources  (explosion  *  tectonic  rv  *.  * 


lig.  4 . b - - Theorct  ical 
sional  arrivals  from 
CIT  109  Low  Q  model, 


spectra  of  tectonically  generated 
the  Handley  event  at  1  distance  of 
explosion  model  T-l,  source  deptli  1 


comp res  - 
4 On (>  Km. 
.2  km . 
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transformation,  the  time  domain  synthetic  signal  pulses  are 
obtained . 

Figure  4.8  shows  individual  teleseismic  compress iona  1 
pulses  with  those  on  the  left  corresponding  to  true  ground 
motion  and  those  on  the  right  the  correspond ing  phases  as 
recorded  on  the  LRSM  short  period  seismometer.  The  seismo¬ 
meter,  which  essentially  measures  particle  velocity,  has  a 
dramatic  effect  on  the  wave  form. 

Phase  Summation 

The  process  of  summation  of  the  individual  phases  (for 
the  vertical  component  of  motion)  in  the  time  domain  is  illus¬ 
trated  in  Fig.  4.9.  The  composite  signal  shown  in  Fig.  4.9f 
corresponds  to  the  first  arrival  mantle  phase  at  4066  km.  It 
is  probable  that  this  entire  wave  form  would  be  taken  to  re¬ 
present  a  single  P  wave  and,  for  example,  the  'first  arrival" 

P  wave  spectrum  would  be  computed  from  it,  giving  the  spectrum 
of  Fig.  4 .  G  cl .  On  the  other  hand,  it  is  clear  that  the  first 
motion  is  appropriate  to  the  direct  explosive  P  wave  from  the 
source  and  the  m^  value  measured  from  peak  to  peak  during 
the  first  cycle  would  be  unaffected  by  pP  reflections  and 
tectonic  release,  even  if  the  tectonic  release  were  several 
times  larger  than  shown  here. 

Figure  4.10  shows  the  expected  time  domain  composite 
signals  with  and  without  tectonic  release  effects.  The  re¬ 
sults  of  Fig.  4.9  and  4.10  are  typical  in  that  tectonic  re¬ 
lease  may  have  an  appreciable  effect  on  the  signal  shape  but, 
in  general,  does  not  affect  m^.  This  is  because  the  major 
tectonic  phase,  sP.j. ,  arrives  too  late  while  the  earlier 
arrivi  ng  P.p  and  pP-j,  are  too  small  to  change  m^ .  Also 
shown  in  Fig.  4.10  are  synthetic  seismograms  including  a 
second  mantle  arrival.  This  illustrates  our  ability  to  gener¬ 
ate  a  long  time  segment  of  the  seismogram  including  all  source 
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(a)  Direct  F  wave  fron  explosion  (P£) 


(b)  Surface  reflected  P  wave  fron  explosion  (pPjr) 


(c)  Direct  P  wave  fron  tectonic  release  (P^) 


(d)  Reflected  P  wave  fron  tectonic  release  (pP^) 


(e)  Surface  converted  S  to  P  wave  fron  tectonic 
release  (sP^) 


1  2  3  4  !>  6 

(f)  Conposite  first  mntle  arrival  series. 


F  ig.  4 . 9 -  - The o re t i ca  1  compressional  (P)  wave  phases  generated  bv 
the  Ilar.dley  explosion  and  associated  tectonic  release.  The 
Handle)  event  flOoO  kT)  modeled  by  the  explosion  source  model  T-.l 
Prestrc-ss  (o15)  for  tectonic  release  taken  to  be  65  bars,  shattci 
zone  radius  Rq  =  750  meters.  Earth  structure  C1T  109- -Low  0 
model.  Distance  4066  km,  azimuth  30°.  Vertical  component  MIS'  I 
short  period  seismometer,  source  depth  1.2  km. 
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effects  as  well  as  free  surface  reflections.  further,  we 
routinely  generate  the  S  (body)  wave  portion  of  the  seismo¬ 
gram  and  surface  waves  (Section  III)  as  well  as  the  P  waves, 
so  that  most  of  the  complex  signal  train  generated  by  an  event 
is  synthesized  in  a  series  of  rather  fast  computations. 

lypical  Seismograms 

lo  ill u st; ate  the  degree  of  agreement  that  is  possible 
with  our  computational  model,  we  compared  the  observed  signal 
from  the  Bilby  event  (  200  kt  in  tuff)  at  LRSM  station  I ! .MM II 
(4000  km  from  the  NTS  site)  to  a  synthetic  seismogram  in  Fig. 
4.  II.*  For  this  computation  the  tuff  model  T-3,  which  has  a 
displacement  pulse  somewhat  larger  than  that  of  T-4  (Fig. 

4.21,  was  used.  While  the  amplitude  scale  of  the  observed 
record  is  not  known,  the  LRSM  reported  A/T  of  107  millimic¬ 
rons  for  IINMF.  'ompares  favorably  with  the  A/T  of  127  milli¬ 
microns  for  the  first  cycle  of  the  synthetic  seismogram. 

I  he  S  wave  portion  of  the  synthetic  seismogram  for 
the  event  ju  t  described  is  portrayed  in  Figs.  4.12  and  4.15. 

In  Fig.  4.12  the  individi:-!  phases  are  shown  along  with  their 
summation,  the  first  mantle  arrival,  for  both  vertical  and 
tangential  motion.  The  fir*-r  two  mantle  arrivals  for  the 
vert i ca .  and  tangential  components  of  motion  are  shown  in 
Fig.  4.13. 


*The  synthetic  seismograms  of  earlier  figures  were  hand  drawn 
* r°m _ computer  plots;  hence  the  smooth  curves.  This  and 
ceeding  synthetic  seismograms  were  drawn  by  an  SC-4020 
Ihe  small  oscillations  at  the  beginning  of  the  record 
introduced  by  the  numerical  Fourier  transform. 


suc- 

plotter, 
a  re 
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(a)  lINMK  LRS'I  SPZ  record  for  Kilby. 


(h)  Vertical  displacement  for  explosion  source  model  1-3  with 
tectonic  release  (prestress  o\2  =  015  =  100  bars,  shatter 
zone  radius  Rq  =  0.6  km)  propagated  with  earth  model  C IT 
109  bow  Q,  source  depth  0.7  km. 


1  ig.  4 . 1 1  - - Compa r i son  of  observed  and  theoretical  seismograms 
for  the  Kilby  event  at  ld.'Mb,  4006  from  ground  zero 


( i.i )  Direct  S-wa vc  iron  tectonic  .source 


10.0  II. 0  II. 0 


(!■>)  Reflected  S-wavc  from  tectonic  s 


source 


1  1  8  •  i.  1  2  -  -  Lon  s  t  met  i  on  of  the  first  m  a  n  1 1  n  nrviini  c  ,  .  r 
vert ica !  displacement  (left  column),  and  tangential  displacement 
tl,  respect  to  y  round  zero  friyl.t  colon,,;)  .  I)  i  s„  1  nee, on  ! 
J'p, laved  n,  null, microns  ami  times  in  seconds-'  he  eve,,?  , 
meters  arc  given  in  fig.  4.11.  1,11,1 
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event  described  in  F  ig'.'Vl  ” 1 *  C;)‘” 1  s  4060  kni  r°r  the  Bilbv 
'}}} d  times  are  in  seconds.  Zero  on  ,  ?°ntS  arC  in  imicrons', 
f,*°  scconds  from  detonation  time  °  SCalc  corrc*Ponds  to 
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The 


r.oct  of  IT,  r  t  li  Structure  on  the  I  c  I  esc  i  sm  i ,-  |ic.col-d 


'  01  rCli'ihl°  tclosoismic  •!  i,l  p  I  i  t  lulc  predictions,  one 

"  1S '  ,avc  tlu’  '"ost  ;lccl,r;'t0  earth  model  aval  ]  able.  The  ;,rri- 
tines  of  the  various  mantle  arrivals  are  entirely  de¬ 
pendent  on  the  model  used.  Further,  since  the  earth' acts  as 
•'  puss  filter,  the  amplitude  and  shape  of  the  arriving 
Piuses  are  stron8ly  dependent  on  the  attenuation  and  disper¬ 
sive  „  rope  rt  i  es  of  the  earth  model  chosen.  Therefore,  the 
earth  model  must,  at  .a  minimum,  specify  the  averaSe  compres- 

.  "ll  ,UU  Slc:lr  eelocit.ics  Kith  depth,  as  Kell  as  the  dis¬ 
sipation  (unction,  q,  as  a  function  of  depth. 


h-ivcrcd  1 : a  r t h  Model 


lor  the  Kork  of  tnis  report  a  layered  earth  model  pro¬ 
posed  hy  Archan, beau,  Minn  and  Lambert'23)  ,vas  used.  This 
model,  known  as  CIT  HIP,  is  appropriate,  in  terms  of  the  ' 
upper  mantle  and  crustal  structure  in  the  upper  2110  km,  to 
10  asm  and  Range  Province  of  the  Western  United  States 
clow  about  200  kn,  lateral  variations  in  structure  witlun 
ie  earth  are  such  that  the  velocity  structure  is  a  reason- 
u>  e  approximation  to  a  mean  earth  velocity  profile  In 

,!**•.  4'U.  1'"d  4JS  the  Volocit>'  0«Pth  profiles  for  compros  - 
I  ona,  and  shear  naves  arc  shown.  The  various  upper  mantle 
•  -scon  nun  ties  including  the  low  velocity  zone  and  especially 
ie  velocity  transitions  at  400  and  650  km  give  rise  to  im¬ 
portant  later  arriving  body  wave  phases  at  telcscismic  dis- 
tunces.  I  he  travel  time  curves  showing  the  timing  of  all 
.•■untie  P  wave  arrivals  for  the  CIT  100  model  are  shown  in 
is-  4.16,  along  with  observations  from  the  Kilby  event.  This 
enure  shows  that  the  first  arrivals  may  often  be  followed, 

""  SLC0lul  or  two>  by  u  second  arrival  and  within  5  or 
seconds  by  as  many  as  3  or  4  arrivals.  In  some  distance 
ranges  these  may  be  as  large  er  larger  than  the  true  first 
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arrival  and  strongly  affect  the  character  of  the  observed 
seismic  recording. 

I  he  dissipation  function  Q  attenuates  and  disperses 
the  telcsc ismic  pulse.  Several  Q  models  were  investigated 
during  the  course  of  this  study.  In  Figs.  4.17  and  4. IS 

several  Q  depth  profiles  for  the  CIT  F  and  CIT  S  wave  models 
arc  shown. 

transfer  Function 

For  the  CIT  109P  0.75  Q  model,  the  amplitude  of  the 
transfer  function,  |  (or)  [  (see  liq .  (4.4))  for  to  =  2  tt  ,  is 
plotted  versus  distance  in  Fig.  4.19.  This  transfer  func¬ 
tion  was  calculated  for  a  source  depth  of  700  meters.  Plots 
°*'  I  •■''p  ( - T’ )  |  for  each  of  the  first  three  arrivals  (at  100 
km  intervals)  are  shown  separately  on  the  figure.  Plots 
for  the  other  Q  models,  CIT  109P  High  Q  and  CIT  109P  how 
Q,  take  a  similar  form,  merely  being  shifted  above  or  be¬ 
low  the  plots  of  Fig,  4.19. 

The  tiansfer  functions  for  other  frequency  components 
of  interest  for  teleseismic  P  waves  exhibit  a  distance 
dependence  similar  to  that  of  the  u  =  2 tt  component,  there¬ 
fore  Fig.  4.19  is  essentially  a  plot  of  amplitude  versus 
distance  for  P  wave  arrivals.  The  travel  time  curves  (Fig. 
4.19)  together  with  the  transfer  function  curves  (Fig.  4.19) 
give  the  relative  amplitudes  and  travel  times  of  the  P  wave 
mantle  arrivals  at  teleseismic  distances.  Two  examples  of 
the  resulting  synthetic  seismograms  are  shown  in  Fig.  4.20. 

•  «ote  that  at  3000  km  the  second  two  arrivals  interfere, 
greatly  complicating  the  record. 
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(a)  First  three  mantle  arrivals  at  2500  km.  Zero  on  the 
time  scale  is  299.2  seconds  from  detonation  time. 
Seismogram  appropriate  to  the  Bilby  event  with  tectonic 
release- -see  Figs.  4.31  and  4.33.) 


(b)  First  three  mantle  arrivals  at  3000  km.  Zero  on  the 
time  scale  is  341.9  seconds  from  detonation  time. 
(Seismogram  appropriate  to  the  Rulison  event  (no  tec¬ 
tonic  release) - -see  Fig.  4.39.) 


Fig.  4. 20- -Two  synthetic  seismograms  illustrating  the  relative 
amplitudes  and  arrival  times  of  the  first  three 
mantle  arrivals.  The  displacements  are  measured  in 
millimicrons. 
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The  1  f feet  of  Q 

Tlie  dissipation  function  Q  attenuates  the  traveling 
T  wave.  Since  our  earth  model  is  dispersive,  0  also  affects 
1  he  shape  of  the  arriving  wave  form,  with  decreasing  Q 
causing  increased  spreading  of  the  arriving  pulse.  This 
effect  is  illustrated  in  Fig.  4.21  where  the  first  mantle 
arrival  from  an  explosion  (with  no  tectonic  release)  is 
depicted  for  the  three  Q  models.  The  quantities  A.,  T.  ,  A_ , 

T-i ,  which  are  defined  in  Fig.  4.21(a),  are  the  amplitudes 
(in  millimicrons)  and  periods  (in  seconds)  of  interest  for 
calculating  m,  . 

The  attenuating  and  dispersing  effects  of  Q  arc  ap¬ 
parent  from  the  three  records  of  the  direct  P  wave.  Since 
the  explosion  depth  is  rather  shallow,  the  direct  and  free 
surface  reflected  P  wave  phases  exhibit  considerable  inter¬ 
ference  with  varying  effects  on  the  seismograph  record,  de¬ 
pending  on  Q. 

Calculation  of  m^ 

The  synthetic  seismograms  of  Fig.  4.21  are  appropri¬ 
ate  to  the  Benioff  short  period  instruments  of  the  LRS'l 
stations.  As  outlined  in  the  LRSM  shot  reports  (e.g.,  Ref. 
[19]),  m^  is  picked  from  the  cycle  of  maximum  amplitude 
among  the  first  few  cycles  of  the  record.  As  will  be  dis¬ 
cussed  in  subsequent  sections,  only  the  first  cycle  is  directly 
related  to  the  explosion  yield  with  later  cycles  being  affected 
by  free  surface  reflected  phases  and,  possibly,  tectonic  re¬ 
lease.  However,  for  comparison  of  our  results  with  observa¬ 
tions,  wre  must  attempt  to  closely  approximate  the  LRS'I  proce¬ 
dure  for  picking  m^.  For  records  like  those  of  Fig.  4.21  it 
is  likely  that  the  second  cycle  (of  amplitude  A-,)  would  be 

L. 

used. 
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(a)  CIT  109P  High  Q.  Aj  =  187,  T1  =  0.6,  ,\9  =  252  ,  T?  =  l, 


1.0  19 


•o  1.0  7.0  3.0  4.0 


(h)  CIT  10!IP  0.75  Q.  Aj  =  168,  Tj  =  0.5,  A,  =  136,  T,  =  1.1 
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(c)  CIT  109P  Low  Q.  A:  =  54,  T1  =  0.75,  A?  =  70,  T?  =  1.1 


Fig.  4.21--Effect  of  Q  on  the  direct  P  phase  (left  column) 
and  first  mantle  arrival  at  4066  km  (sum  of  direct  and  free 
surface  reflected  P  waves).  Source  T-4,  DOB  700  m,  yield 
200  kt.  Displacements  are  in  millimicrons. 
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The  body  wave  magnitude  is  calculated  from 

m^  =  log  A^/'IT  +  B,  (4.5) 


is  the  distance  factor  3.54  for  4066  kmj  and 


-ig' 


4.21 

(2) 


The  values  of  m,  for  the 


)  and  second  (m,  ^  J  )  cycles  are  listed  in  Table 
D  ( 1 1 

v  J  for  each  0  model  is  due 


m, 


where  1 

V,  i' .  are  defined  on 
first 

4.1.  The  higher  value  of 

U 

to  the  smaller  period  for  this  cycle.  The  source  parameters 
for  Fig.  4.2]  are  appropriate  to  the  Bilby  event,  which  had 
an  LRS.M  assigned  m^  of  5.8,  placing  the  values  of  Table  4,1. 
well  within  the  acceptable  range. 


TABLE  4.1 

COMPARISON  OF  mb  FOR  RECORDS  OF  FIGURE  4.20 


Q  Model 

mb(1) 

(2) 

mb 

High  q 

6.03 

5.91 

0.7  5  q 

6.07 

5.79 

Low  q 

5.40 

5,34 

Summarv  of  Ea:*th  Structure  Effects 

-  r  .  - -  .  — —  i  ■■  


The  accuracy  of  our  predictions  of  and  the  form  of 
the  teleseismic  record  are  quite  strongly  affected  by  the 
velocity  and  Q  structures  used.  The  CIT  109  P  and  S  velocity 
models  appear  to  give  good  results,  but  improved  detailed 
agreement  with  observations  could  be  obtained  by  refinement 
or  replacement  of  these  models.  For  example,  it  would  be 
straightforward  to  include  two  layered  earth  models  in  our 
computations;  one  for  propagating  the  ray  from  the  source 
region  to  its  deepest  penetration,  and  one  for  propagating 
the  ray  from  its  turning  point  to  the  receiver,  which  may 
be  located  on  a  quite  different  earth  structure.  Further, 
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since  crustal  reverberations  have  a  significant  effect,  these 
also  could  be  modeled  for  receivers  of  interest. 

4.3.3  The  hffect  of  Depth  of  Burial  on  the  Teleseismic  Record 

The  burial  depth  of  the  explosion  source  affects  the 
teleseismic  record  in  several  fundamental  ways.  By  far  the 
most  important  of  these  is  the  interference  between  the  direct 
(P)  and  free  surface  reflected  (pP)  waves  leaving  the  source 
vicinity.  Other  depth  of  burial  (DOB)  effects  include  the 
variation  of  local  material  properties  and  overburden  pres¬ 
sure.  Finally,  great  DOB  changes  could  significantly  alter 
the  transfer  function  for  the  source  to  receiver  travel  path. 

Material  property  changes  are  properly  considered  a 
source  effect  since  the  source  spectrum  is  directly  a  function 
of  the  material  properties  in  the  source  region.  The  effect 
of  important  source  material  parameters  on  the  source  spect¬ 
rum  is  discussed  in  Section  II.  Further,  it  is  to  be  ex¬ 
pected  that  the  average  source- to- surface  P  and  S  wave  velo¬ 
cities  may  also  change  with  depth,  thereby  changing  the 
apparent  depth  of  burial.  For  example,  a  10  percent  decrease 
in  the  average  P  wave  velocity  results  in  an  11  percent  in¬ 
crease  in  the  arrival  time  lag  between  P  and  pP,  having  the 
effect  of  increased  depth  of  burial. 

The  overburden  pressure  P^  increases  with  increasing 
depth,  having  the  apparent  effect  of  slightly  decreasing  the 
vie  Id .  In  Fig.  4.22  the  amplitudes  of  the  RYP  transforms 
l'i'(to)  ,  from  three  SKIPPER  calculations  are  compared.  The 
three  source  models,  T-4,  T-5  and  T-6  in  Table  2.2,  Section 
II,  are  identical  except  for  overburden  pressure.  Kith  in¬ 
creasing  depth  of  burial  we  see  that  the  RVP  transform  ex¬ 
hibits  the  behavior  predicted  by  the  scaling  laws  (Table  2.1) 
for  decreasing  yield;  that  is,  reduced  ip ( 0 )  and  a  spectral 
shift  to  higher  frequencies.  These  effects  are  quite  minor, 


98 


R-2034 


1.75-10 


1.5*10 


1 .25x10 


7.5*10 


2.5x10 


T--1 

I 

1 

• 

30 

Frequency  (lie) 


Fig.  4 . 22--Reduced  velocity  potential  transform  amplitude 
versus  frequency  for  tuff  source  functions  T-4,  T-5  and 
T-6.  Scale  is  appropriate  to  a  yield  of  D.02  kt. 
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however,  as  a  more  than  twofold  increase  in  depth  results  in 
an  apparent  yield  decrease  of  about  15  percent  (based  on 
>‘(0)). 

Interference  Between  P  and  pP 

The  most  striking  feature  of  the  body  wave  spectra  from 
an  explosion  is  the  peaking  and  scalloping  of  the  spectrum  due 
to  interference  between  the  direct  and  reflected  waves, 
examples  were  shown  in  Figs.  4.5  and  4.6.  The  effect  of  the 
free  surface  on  the  spectrum  of  the  observed  first  arrival  P 
wave  train  can  be  demonstrated  as  follows: 

If 

P  =  F  Ct ) ,  (4.6a) 

then 

pP  =  -AF(t  -  Tp,  (4.6b) 

where 


A  free  surface  reflection  coefficient  (f  l)  , 

T  „  2!!/v  , 
o  -  p* 

II  =  depth  of  burial  (DOB), 

Vp  =  compressional  wave  velocity. 

In  the  frequency  domain 

a  a  /  "iwT  \ 

P  +  pP  =  F(u>)  \  1  -  Ae  °/, 
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and 


|P  +  pP 


FO) 


1  -  2 A  cos  (gjT  ) 


1/2 


[4.7) 


(u>) 


1  -  A 


1+A 


U)  - 


a)  - 


♦ 

-  To 


.  (2n  + 1 )  tt 

_  — j - »  n 


(4.3) 


=  012 
v'  >  1  |  •  •  • 


Equation  (4.8)  shows  the  modulation  of  the  amplitude 
spectrum  by  the  free  surface  reflection,  with  the  spectral 
troughs  and  peaks  occurring  (alternately)  at  frequencies  of 


f  -  ,  n  =  0,1,2,  ...  .  (4.9) 

The  effect  of  Eq.  (4.9)  is  shown  in  Fig.  4.23  where 

/V  /V  /\  /\ 

I P I ,  j  pP |  and  ]P  +  pP |  are  plotted  for  two  hypothetical  yields 
and  burial  depths.  The  two  calculations  were  run  using  the 
same  source  function  and  therefore  the  same  overburden  pres¬ 
sure  (appropriate  to  a  DOB  of  700  m) .  However,  using  the 
proper  overburden  pressure  for  the  deeper  depth  would  result 
in  only  a  slight  downward  shift  of  the  three  curves  for  1400 
kt.  The  material  properties  in  the  source  region  were  also 
held  constant  for  both  calculations,  with  the  source-to-sur- 
face  compre s s ional  wave  velocity  being  2.4  km/sec. 

To  demonstrate  that  the  spectral  shift  of  Fig.  4.23  is 

almost  entirely  a  DOB  effect,  in  Fig.  4.24  the  vertical  com- 

/*\  /*>. 

ponent  of  |P  +  pP|  is  plotted  for  six  yields  ranging  from  1 
to  800  kt  with  the  source- to- surface  travel  time  remaining 
constant . 
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Vertical  Displacement  (nu/cps) 


-Vertical  component  P  wave  spec 
distance  of  3000  km  for  RDP  s 
yields  and  burial  depths. 


102 


R  -  2  0  5  4 


Tlie  effect  of  free  surface  reflections  on  the  tele- 
seismic  record  is  best  examined  by  fixing  all  parameters 
(yield,  source  spectrum,  source- to- surface  velocity)  but  the 
depth  of  burial  and  studying  synthetic  seismograms  due  to 
the  explosion  only  (no  tectonic  release).  Since  the  trans¬ 
fer  function  for  the  source-receiver  ray  path  has  equal 
effect  on  P  and  pP,  only  the  arrival  time  lag  between  the 
two  is  significant. 


In  Fig.  4.25,  synthetic  seismograms  for  depths  of  700 
and  1200  meters  are  compared.*  Using  Eq.  (4.9),  the  first 
spectral  window  for  the  events  of  4.21a  and  4.21b  occurs  at 
1.7  and  1.0  Hz  respectively.  Using  the  same  source  para¬ 
meters  but  a  different  earth  model,  the  depths  of  700  and 
2573  meters  are  compared  in  Fig.  4.26.  At  the  greater  depth 
pP  is  so  delayed  as  to  appear  as  a  later  mantle  arrival. 

The  calculation  of  from  records  such  as  those  of 
Figs.  4.25  and  4.26  is  ambiguous.  As  briefly  discussed  in 
Section  4.3.2,  mb  in  calculated  from  the  cycle  of  maximum 
amplitude  from  among  the  first  few  cycles  of  the  record. 

The  formula  (4.5)  is  then  applied.  Using  this  LRSM  procedure 
mb  would  almost  certainly  be  calculated  from  the  second  cycle 
on  the  shallow  shots  (Figs.  4.25b  and  4.26b)  and  from  the 
first  cycle  on  the  deeper  events.  The  mb  values  for  these 
records  are  summarized  in  Table  4.2.  It  should  be  pointed 
out  that  the  mb ^  values  are,  in  a  sense,  anomalous  since 
the  periods  ^  are  cn  the  order  of  0.5  sec,  values  seldom 
observed  for  mb  measurements  at  teleseismic  distances. 


*  See  Figure  4.21  for  definitions  of  A.  and  T. 
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(a)  DOB  =  1200  meters.  =  215,  =  0.5. 
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(b)  DOB  =  700  meters.  =  187,  =  0.6,  A 2  =  232,  T2  -  1. 


Fig.  4 . 25- -Comparison  of  first  arrival  P  waves  (P+pP) ,  as 
recorded  at  4066  km,  for  different  depths  of  burial.  Source 
T- 4 ,  CIT  109P  High  Q  earth  model,  vp  =  2.4  km/sec.  Displace¬ 
ments  are  in  millimicrons. 
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(a)  DCB  -  2573  meters.  Ax  =  230,  T  =  0.5 


(b)  DOB  =  700  meters.  Ax  =  160,  ^  =  0.5,  A 


1?5,  T2 


Fig.  4 -26- -Comparison  of  first  mantle  arrival 
Scarce  T-4,  CIT  109P  0.75  Q  earth  model 
Placements  are  measured  in  millimicrons!  p 


waves  at 
km/sec . 
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TABLE  4.2 

COMPARISON  OF  mb  FOR  RECORDS  OF  FIGURES  4.24  AND  4.25 


R- 2034 


Summary  of  Depth  of  Burial  Effects 

Although  DOB  has  some  minor  influence  on  the  source 
spectrum,  the  dominant  DOB  effect  is  the  shifting  of  the 
arrival  time  of  pP  with  respect  to  P.  For  shallow  depths 
the  two  phases  may  combine  to  significantly  increase  the 
maximum  amplitude  on  the  record.  This  maximum  amplitude, 
from  which  mfe  measurements  are  made,  is  therefore  dependent 
not  only  on  explosion  yield  but  also  on  DOB.  Below  a  certain 
depth  (about  1000  m  for  the  tuff  example  treated  here),  the 
two  pulses  disengage  on  the  seismogram  and  the  maximum 
amplitude  is  then  directly  dependent  on  yield. 

4  *  3  * 4  .The  Effect  of  Yield  on  the  Teleseismic  Record 

Taken  by  itself,  the  effect  of  the  explosion  yield  is 
essentially  restricted  to  scaling  the  record  amplitudes 
directly  with  yield.  In  practice  a  change  in  yield  is 
generally  associated  with  a  change  in  the  depth  of  burial 
and,  possibly,  other  parameters  of  the  event,  thus  obscuring 
this  elementary  scaling  relationship. 

In  Fig.  4.24  the  total  displacement  spectra  due  to  a 
series  of  explosions  over  a  wide  range  of  yields  were  dis¬ 
played.  For  these  events  the  source  environment  and  DOB 
were  held  constant;  therefore,  the  effect  of  yield  is  iso¬ 
lated  and  is  restricted,  aside  from  the  amplitude  scaling, 
to  a  very  slight  shift  of  the  spectrum.  This  result  is  a 
s traight forward  consequence  of  the  properties  of  the  earth, 
as  will  be  explained  below. 

Source  Spectrum  Changes  with  Yield 

The  calculation  of  the  source  function  by  SKIPPFR  is 
carried  out  at  a  nominal  yield  of  0.02  kt,  then  scaled  to 
the  yield  of  interest.  The  far  field  displacement  spectrum 
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is  directly  proportional  to  the  reduced  velocity  potential 
CRVP)  transform,  Hu)  (see  Eq.  (4.1)).  Therefore,  the 
scaling  laws  of  interest  are  that  ^(ui)  a,  and  w  v 
where  li  is  the  explosion  yield. 

In  Fig.  4.27  |i^(ui)|  for  the  SKIPPER  calculation  T-4 

(tuff)  is  shown.  The  spectral  shape  is  typical  of  our  source 
calculations.  The  amplitudes  and  frequencies  given  on  the 
axes  are  for  a  yield  of  0.02  kt.  Since  the  amplitude  scales 
directly  with  yield,  appropriate  amplitude  values  for  any 
desired  yield  are  quickly  obtained.  Of  greater  interest  is 
the  observation  that  the  spectral  shape  changes  with  yield. 

Ihis  is  indicated  by  the  additional  frequency  scales  on  the 
figure  for  yields  of  200,  500  and  1000  kt .  As  yield  increases, 
the  1  1 1  z  point  on  the  spectrum  moves  down  the  flat  part  of  the 
spectrum,  through  the  trough,  etc.  Therefore,  the  displace¬ 
ment  spectrum  in  the  source  region  varies  significantly  with 
yield . 

Propagation  of  the  source  spectrum  to  teleseismic 
distances  is  accomplished  by  multiplying  by  the  transfer 
function  Eq.  (4.4).  The  earth  acts  as  a  low  pass  filter  with 
the  filter  attenuation  given  by  Eq.  (4.2).  In  Fig.  4.28  the 
attenuation  function  exp  [-u)T/2Qp]  is  plotted  for  several 
typical  ray  paths.  From  this  figure  it  is  seen  that  the 
attenuation  is  quite  strong,  so  strong  that  details  of  the 
source  spectrum  are  essentially  washed  out.  That  is,  even 
such  differing  source  -pectra  as  those  indicated  in  Fig. 

4.2/  result  in  quite  similar  spectra  at  teleseismic  distances. 
The  latter  assertion  is  born  out  by  Fig.  4.29  where  direct  P 
wave  spectra  at  4066  km  are  compared  for  several  yields.  The 
spectrum  for  0.02  kt  has  been  scaled  by  104  to  juxtapose  it 
with  the  others. 


109 


At  tcniin 


<000  kn 


2334  k- 


Frequency  (I!:) 


Fig.  4.28- -The  attenuation  function  exp [ -wT/20n]  for  three  rn\ 
paths  in  the  model  CIT  109  0.75  0.  All  other  attenuation 
functions  for  this  model  fall  in  the  shaded  area. 


Vertical  Displacement 


Frequency  (l!z) 


Fig.  4 . 29- -Spectral  amplitude  of  the  direct  P  wave  at  4066  km 

nnAV16  source  spectra  of  Fig.  4.27.  The  amplitudes  for 
0.02  kT  have  been  multiplied  by  104. 
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Scaling  of  nh  l.’ith  Yield,  V,' 

In  \  ieu  of  the  above  discussion,  m^  versus  log  If  curves 
>ue  expected  to  have  slope  ~  1.  Fixing  all  other  parameters, 
the  teleseism ic  records  of  the  first  mantle  arrival  due  to 
three  events  of  different  yields  are  shown  in  Fig.  4.50.  \'o 

differences  in  the  pulse  shapes  for  the  three  yields  are 
appa rent . 

]•  igure  4.30d  is  essentially  a  plot  of  m^  versus  log 
W.  If,  following  Basham  and  Horner, we  uere  to  define 

mb  =  n’h  +  log  T»  (4.10) 

and  plot  m£  versus  log  h ,  we  would  obtain  a  parallel  line. 
Although  the  source  spectrum  is  quite  different  for  the  three 
yields  (Fig.  4.27),  only  a  slight  deviation  from  unit  slope 
is  apparent  in  the  plot.  A  linear  fit  to  the  data  would  have 
a  slope  of  x  0.975. 

4.3.5  Tectonic  Stress  Release 

Tectonic  stress  release  is  recognized  as  the  most  likely 
cause  of  the  anomalous  part  of  the  radiation  field  from  large 
underground  explosions  (e.g.,  Ref.  [13-16]).  The  theory  0f 
stiess  wave  radiation  from  explosions  in  prestressed  media 
formulated  by  Archambeau ^ ^ ^  is  used  to  compute  the  tectonic 
release  contribution  to  the  radiation  field.  Typical  cal¬ 
culations  including  tectonic  release  effects  were  presented 
in  Section  4.3.1,  where  it  was  observed  that,  in  general, 
tectonic  release  does  not  affect  m^ ,  but  may  have  an  appreci¬ 
able  effect  on  the  observed  wave  form. 

The  tectonic  release  component  of  the  source  is  speci¬ 
fied  in  terms  of  a  quadrupole  point  source.  The  important 
parameters  controlling  this  source  are  the  magnitude  and 
orientation  of  the  prestress  and  the  radius  (R  )  of  the 
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shatter  zone  created  by  the  explosion  shock  wave.  The  dis¬ 
placement  spectral  amplitudes  scale  directly  with  the  pre¬ 
stress  and  exhibit  the  azimuthal  dependence  characteristic 
of  a  quadrupole  source.  The  amplitudes  also  scale  with  the 
cube  of  the  shutter  zone  radius.  Since  R  is  expected  to 
scale  as  IV  ,  the  latter  is  equivalent  tc  tectonic  release 
amplitude  scaling  with  yield. 

Tectonic  Release  Effects  on  SPZ  Records 

The  theory  of  tectonic  stress  release  from  explosions 
was  developed  to  account  for  observations  of  anomalous  shear 
wave  radiation  from  large  underground  explosions.  Most 
studies  of  this  phenomenon  are  primarily  concerned  with  sur¬ 
face  waves,  however,  let  us  examine  the  possible  contribution 
of  tectonic  release  to  measurements  of  shore  period  vertical 
motion  (SPZ).  In  Fig.  4.31  the  true  ground  motion  for  the 
phases  which  make  up  the  first  mantle  arrival  from  an  explo¬ 
sion  with  tectonic  release  are  shown,  as  well  as  the  records 
of  this  motion  from  a  Benioff  short  period  instrument  (see 
Fig.  4.7  for  the  instrument  transfer  function).  The  records 
are  appropriate  for  an  epicentral  distance  of  4066  km  and  an 
azimuth  of  45  degrees.  The  conventions  for  defining  pre¬ 
stress  and  source- to-receiver  azimuth  are  shown  in  Fig.  4.32. 
From  Fig.  4.31,  note  that  the  first  P-wave  motion  from  tec¬ 
tonic  release  is  rarefaction,  as  expected.  From  symmetry, 
the  records  of  Fig.  4.31  also  apply  for  a  prestress  of 
al2  ~  "160  bars  and  an  azimuth  of  225  degrees. 

Ihe  effect  of  the  LRSM  SP  instrument  is  clearly  indic¬ 
ated  in  Fig.  4.31.  It  is  interesting  to  note  that  for  the 
simple  direct  P  wave  pulse  from  the  explosion,  the  peak-to- 
peak  amplitude  on  the  instrument  record  is  nearly  the  same 
as  the  peak  true  ground  motion.  The  more  complicated  tec¬ 
tonic  release  wave  forms  do  not  exhibit  this  relationship, 
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(a)  Direct  P  wav.1  from  the  explosion. 


(b)  Reflected  P  wave  from  the  explosion. 


Fig.  4.31--P  wave  phases  at  a  distance  of  2500  km  and  an 
azimuth  of  45  degrees  from  a  source  consisting  of  an  explo¬ 
sion  plus  tectonic  release.  Given  are  the  true  ground  motion 
(right  column)  and  seismometer  output  from  LRSM  SP  instrument. 
(Explosn  i  source  T-4,  yield  =  200  kt,  DOB  =  700  m,  Ro  =  600  m, 
prestress  =  160  bars,  earth  structure  CIT  109P  0.75Q.) 

The  displacements  are  indicated  in  millimicrons. 
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x-  (vertical) 


Fig.  4.32  --Coordinate  system  for  defining  the  prestress  (arrows) 
and  source- to-receiver  azimuth  (9).  The  first  motion  in  the 
four  quadrants  due  to  an  explosion  in  a  medium  with  prestress 
ax2  >  0  is  also  indicated.  Definitions  of  013  and  023  are 
consistant  with  that  of  o^* 
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clue  mainly  to  the  instrument  filtering  of  the  long  period 
components  of  these  pulses.  The  long  period  lobes  on  the 
tectonic  release  ground  motion  records  (caused  by  finite 
Rs>  see  Fig-  4.4  and  Section  4.3.1)  have  little  effect  on 
the  seismograph  records. 

For  the  event  of  Fig.  4.31,  the  tectonic  P  wave  is 
of  the  same  order  as  the  explosion  P  wave  and  the  tectonic 
sP^  is  the  largest  phase  of  all.  This  is  due  to  the  rather 
large  prestress  and  crush  zone  radius  selected,  as  well  as 
the  station  orientation  with  respect  to  the  prestress  (radia¬ 
tion  pattern  effect).  The  effect  of  this  large  tectonic 
release  component  is  ascertained  by  comparing  first  mantle 
arrival  records  with  and  without  tectonic  release.  In  Fig. 
4.33  the  total  first  arrival  record  (from  the  LRSM  SP 
instrument)  including  only  explosion  effects  is  shown  as 
well  as  the  total  record  with  both  explosion  and  tectonic 
release  phases  included.  In  Table  4.3  the  m^  values  for  the 
two  records  are  compared  (the  distance  correction  for  2500 
km  is  taken  to  be  3.25). 

TABLE  4.3 

COMPARISON  OF  mb  FOR  THE  RECORDS  OF  FIGURE  4.33 


Seismogram 

b 

(2) 

mb 

4.33a  (no  tectonic 

5.62 

5.45 

release ) 

4.33b  (with  tec- 

5.69 

5.55 

tonic  release) 

Clearly,  tectonic  release  has  little  effect  on  mu  for 

b 

this  example.  However,  comparing  the  two  records,  tectonic 
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(a)  First  mantle  arrival  with  no  tectonic  release:  sum  of 
phases  4.31(a)-(b).  A,  =  139,  T,  =  0.6,  =  167, 

x2  =  1.05.  1  1  2 


time  <scci 


(b)  First  mantle  arrival  including  tectonic  release;  sum 
of  phases  4.31(a)-(e).  A,  =  165,  T.  =  0.6,  An  =  199, 

T  =  1  -L  1 


Fig.  4 . 33- -Comparison  of  P  wave  first  mantle  arrival  records 
for  the  event  of  Fig.  4. 31, with  and  without  tectonic  release. 
The  first  three  mantle  arrivals  as  recorded  at  2500  km  are 
shown  in  part  (c).  All  displacements  are  indicated  in 
mill imicrons . 
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release  does  add  considerable  complexity,  especially  via  the 
large  and  late  arriving  s?T  phase. 

In  Fig.  4.33c  the  complete  record  of  the  first  three 
P  wave  mantle  arrivals  is  shown.  At  this  distance,  2500  km, 
the  first  mantle  arrival  is  unusually  small  compared  to  later 
arrivals.  This  suggests  another  source  of  difficulty  in 
choosing  m^ ,  especially  for  low  yield  events,  in  that  a  later 
arriving  phase  might  be  selected  for  computation  of  m^ . 

4 • 4  THE  DEPENDENCE  OF  mb  ON  EARTH  STRUCTURE  AND  HIE 
SOURCE  ENVIRONMENT 

4.4.1  General  Remarks 

In  previous  sections  we  have  seen  how  the  various  para¬ 
meters  characterizing  an  explosion  event  influence  teleseismic 
amplitude  measurements.  These  parameters  include  material 
properties  in  the  source  vicinity  (determining  the  explosion 
and  tectonic  release  spectra),  the  depth  of  burial,  the  yield, 
and  the  velocity  and  attenuation  (Q)  properties  of  the  earth 
structure  between  source  and  receiver.  Using  this  information, 
we  proceed  to  study  a  series  of  events  in  tuff,  granite  and 
sandstone  media. 

4 • 4 • 2  Underground  Explosions  in  NTS  Tuff 

All  calculations  presented  in  Section  4.3  pertained  to 
NTS  tuff.  Reference  was  previously  made  to  such  tuff  events 
as  Bilby  and  Handley,  with  which  our  computational  procedure 
gave  excellent  agreement.  It  has  been  postulated  that  signi¬ 
ficant  amounts  of  tectonic  stress  release  accompanied  many 
of  the  tuff  events.  However,  it  was  learned  in  Section  4.3.5 
that  tectonic  release  has  little  effect  on  the  first  two 

cycles  of  the  teleseismic  record,  the  region  in  which  our  m, 

b 

measurements  are  made.  Therefore,  for  simplicity  let  us 
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assume  no  tectonic  release  component  and  construct  synthetic 
seismograms  for  several  events  in  NTS  tufi.. 

Bilby 

This  much  studied  event  was  the  source  of  the  data  on 

which  the  earth  structure  model,  CIT  109,  used  in  this  report 

was  constructed .  ^ The  event  was  a  200  kt  explosion  in 

NTS  tuff  at  a  depth  of  burial  of  700  m.  Synthetic  seismograms 

applicable  to  Bilby  have  been  presented  in  Figs.  4.11,  4.12, 

4.13,  4.21,  4.25b,  4.26b,  and  4.33.  The  results  from 

these  figures  are  summarized  in  Table  4.4.  The  LRSM  assigned 

m,  for  Bilby  was  5.8;  hence,  the  values  of  the  table  show 
b 

good  agreement  with  the  observed. 

Before  directing  our  attention  to  other  tuff  events, 
we  calculate  synthetic  seismograms  for  the  Bilby  event  at  a 
series  of  epicentral  distances.  These  are  shown  in  Fig. 

4.34.  In  Fig.  4.35  the  maximum  amplitudes,  A ^  and  A,  of 
Fig.  4.34,  are  plotted  versus  epicentral  distance.  Also 
shown  are  the  amplitudes  of  the  "principal  phases"  as  re¬ 
corded  by  the  LRSM  network. These  observed  amplitudes 
should  be  the  maximum  amplitudes  recorded  during  the  first 
•few  cycles  of  motion.  Considering  the  wide  scatter  of  this 
data,  our  agreement  seems  quite  satisfactory. 

The  m^  values  for  the  synthetic  seismograms  of  Fig. 
4.34  are  summarized  in  Table  4.5.  Once  again,  our  calcula¬ 
ted  values  show  excellent  agreement  with  the  observed  mb  of 
5.8.  Recall  that  m^  is  that  which  likely  corresponds  to 
the  LRSM  value,  while  m^  reflects  most  accurately  the  ex¬ 
plosion  yield  (relatively  unaffected  by  DOB)  and  is  given 
for  reference. 

Finally,  note  that  the  variation  of  mb  with  distance 
in  Table  4.5  simply  reflects  an  incompatibility  of  our  earth 
structure  model  with  the  LRSM  distance  correction  factors. 
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Fig.  4. 34- -Synthetic  seismograms  for  Bilby  at  a  series  of 
epicentral  distances.  Source  13B,  Yield  200  kt,  BOB  700 
meters,  earth  structure  CIT  109P  0.75  Q. 
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Fig.  4 . 35- -Maximum  amplitude 
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Improvements  to  our  mode1,  and/or  the  distance  factor  cable 
are  necessary  to  remove  this  incompatah il ity . 

Knickerbocker,  Greeley.  Denham 

Synthetic  seismograms  for  three  MTS  tuff  events  of 
widely  differing  yields  are  shown  in  Fig.  4.36.  These 
events,  Knickerbocker,  Greeley  and  Benham,  were  detonated 
at  virtually  the  same  location  at  depths  of  burial  depen¬ 
dent  on  yield.  ^8’24-*  The  synthetic  seismograms  were  con¬ 
structed  for  an  epicentral  distance  of  3000  km  using  the 
earth  structure  model  CIT  109P  0.75  Q  and  the  SkIPPKR 
Source  T-4. 

The  Knickerbocker  event,  like  Bilby  detonated  at  a 
depth  sufficiently  shallow  to  allow  interference  between  P 
and  pP,  gives  a  seismogram  for  which  m^2^  seems  most  com¬ 
parable  to  LRSM  values.  However,  for  the  two  larger  events, 

P  and  pP  arrive  separately  and  no  large  second  cycle  appears 
on  the  records. 

The  m,  values  for  these  three  events  are  summarized  in 
b  .  . 

Table  4.6.  Also  shown  is  m£  defined  from  the  maximum  ampli¬ 
tude  on  the  record  (4. 10),  along  with  average  m£  values  from 
the  Canadian  Seismograph  Network  given  by  Basham  and  Horner. 
Once  again,  our  computed  m^  values  agree  quite  well  with  the 
observed . 

In  Fig.  4.37  the  teleseismic  data  for  the  four  tuff 
events  studied  in  this  section  is  plotted  versus  yield.  The 
slope  of  the  lines  connecting  adjacent  events  is  also  indic¬ 
ated  on  the  figure  and  should  be  compared  with  the  analogous 
data  (Fig.  4.30)  with  depth  of  burial  effects  removed.  Com¬ 
paring  plots  4.37(a)  and  (b) ,  we  see  a  clear  indication  of 
the  sensitivity  of  mb  to  the  measurement  period.  This  period 
is  quite  difficult  to  measure  accurately  on  our  synthetic 
seismogram  and  is  expected  to  introduce  a  random  error. 
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(a)  Knickerbocker  (71  kt  at  630  m)  .  A,  =  43,  1,  =  0.45,  A_  =  56, 
T,  =  1.0.  Z 


■*  1.0  f.O  1,0  <.0 


(b)  Greeley  (825  kt  at  1214  m) .  Ax  =  560,  ^  =  0.6, 


Fig.  4.36--First  mantle  arrival  synthetic  seismograms  for  the 
Knickerbocker,  Greeley  and  Benham  events  in  NTS  tuff.  Epi- 
central  distance  3000  km,  Source  T-4,  earth  structure  model 
CIT  109P  0.75  Q. 


f 


f  mb  and  maximum  amplitude  versus  yield  for  four  ex 
The  slope  of  the  linear  segment  connecting  variou 
cated . 


For  shallow  explosions  the  effect  of  small  changes  in 
DOB  is  significant  because  it  affects  the  measurement  period. 
This  accounts  for  the  seemingly  anomalous  slope  of  0.79  on 

4.37a.  For  deeper  events,  DOB  becomes  insignificant  for  m. 
measurements.  ° 

As  a  final  observation,  mb  data  from  a  series  of 
different  yields,  and  thus  different  burial  depths,  may  in¬ 
clude  data  obtained  from  different  cycles  for  different 
events.  Thus,  caution  must  be  exercised  when  plotting  m, 
versus  log  yield.  For  example,  using  the  maximum  cycle  on 
the  records  of  Fig.  4.36  (as  suggested  for  the  LRSM  network) 
results  in  an  anomalous  bend  in  the  m,  versus  log  IV  or  m* 

versus  log  If  curve  (Fig.  4.37).  b  b 

4,4,3  Underground  Explosions  in  Sandstone 

Several  nuclear  explosions  in  sandstone  have  been  de¬ 
tonated  in  western  Colorado.  In  this  section,  computations 

applicable  to  two  of  these  events,  Rulison  and  Rio  Blanco, 
will  be  presented. 

Rulison 


Rulison  was  a  gas  stimulation  explosion  of  40  kt  de¬ 
tonated  at  a  burial  depth  of  2573  m.  The  source  function 
calculation  for  Rulison  was  described  in  Section  II.  The 
amplitude  of  the  RVP  transform  from  this  calculation  is 
shown  m  Fig.  4.38.  It  is  seen  that  the  source  spectrum  is 

nearly  flat  in  the  region  (<  3  Hz)  which  influences  the  tele- 
seismic  signal. 


In  Fig . 
Rulison  at  an 
The  record  of 
a  rather  large 
to  illustrate 


4.39  synthetic  seismograms  constructed  for 
epicentral  distance  of  2500  km  are  presented. 
Fig.  4.39b  is  from  a  calculation  including 
tectonic  release  component  and  is  included 
the  possible  effect  of  tectonic  release  on  a 
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(a)  First  mantle  arrival  (no  tectonic  release).  A,  =  30, 
T1  =  0.6.  1 


Tint  utei 

(b)  First  mantle  arrival  including  tectonic  release  (shatter 
zone  radius  Rq  =  400  m,  prestress  012  =  256  bars).  A,  = 
46,  Tx  =  0.7,  A3  =  74,  T.  =  1.1.  “  1 


Fig.  4.39--First  mantle  arrival  for  the  Rulison  explosion  with 
and  without  tectonic  release.  Epicentral  distance  2500  km, 
SKIPPER  Source  S-2,  earth  structure  model  CIT  109P  0.75  Q, 
source  -  to- surface  velocities  assumed  identical  to  velocities  in 
the  source  region. 
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very  deep  event.  The  cycle  of  maximum  amplitude  for  this 
record  is  the  third. 

The  computed  m^  and  m£  values  for  Rulison  are  compared 

to  average  observed  values  given  by  Lambert  and  Abner and 

f  1 

Basham  and  Horner L“  J  in  Table  4.7. 

The  agreement  of  computed  and  observed  m^ ;  that  is, 
the  maximum  amplitude  appearing  on  the  record,  is  excellent 
for  the  explosion  without  the  tectonic  component.  The  poorer 
agreement  of  m^  values  then  seems  to  be  due  to  the  relatively 
short  periods  of  the  cycles  from  which  our  measurements 
are  made.  This  is  believed  to  be  a  consequence  of  our  earth 
structure  model,  but  indicates  an  area  which  warrants  further 
investigation. 

The  tectonic  release  component  in  Fig.  4.39b  seems  too 
high  for  the  Rulison  explosion,  but  represents  a  not  im¬ 
plausible  event  for  a  tectonically  stressed  sandstone  region. 

Rio  Blanco 

For  Rio  Blanco  three  30  kt  devices  were  simultaneously 
detonated  in  a  vertical  hole  at  depths  of  1780,  1900  and 
2040  meters.  The  SKIPPER  source  S-l  was  calculated  for  a 
depth  intermediate  to  the  three,  giving  the  |  ijj |  shown  in 
Fig.  4.40.  The  amplitude  scales  on  this  figure  are  appropri¬ 
ate  to  a  yield  of  30  kt  while  frequency  scales  for  both  30 
and  90  kt  are  shown. 

Ignoring  depth  effects,  which  are  expected  to  be  small, 
two  ways  of  treating  this  event  are  suggested.  First,  we 
could  compute  the  teleseismic  signal  due  to  a  30  kt  explosion 
and  multiply  the  far  field  signal  by  three.  Second,  we  could 
treat  the  event  as  a  single  90  kt  explosion.  Examining  Fig. 
4.40,  it  is  apparent  that  the  two  approaches  are  equivalent 
due  to  the  flat  character  of  the  spectrum. 
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Ihc  first  mantle  arrival  portion  of  a  synthetic  seis¬ 
mogram  tor  Rio  Blanco  is  shown  in  fig.  4.41.  The  cycle  of 
maximum  amplitude  on  the  record  is  the  first  which  is  in¬ 
fluenced  only  the  the  direct  P  wave  from  the  explosion.  For 
this  computation  the  body  wave  magnitudes  are 


(4.11) 

mg  =  5.29. 

C 

4-4-4  Underground  Fxplosions  in  Granite 

This  section  includes  an  investigation  of  the  tele- 
seismic  signature  of  explosions  in  granite.  One  specific 
event,  Piledriver,  is  modeled  along  with  a  study  of  magni¬ 
tude  versus  yield  and  tectonic  release  effects. 

Source  Functions 

The  calculation  of  source  functions  for  explosions  in 
granite  is  discussed  in  Section  II,  where  two  SKIPPER  cal¬ 
culations  giving  markedly  different  results  were  presented. 

The  amplitudes  of  the  RVP  transforms  for  the  two  source  func¬ 
tions  are  shown  in  Fig.  4.42. 

Piledriver 


The  Piledriver  event  was  a  61  kt  explosion  at  a  depth 
of  463  meters  in  NTS  granite.  Synthetic  seismograms  using 
both  source  functions  G-l  and  G-2  are  presented  in  Fig.  4,45. 
These  seismograms  show  the  first  mantle  arrival  at  an  epi- 
central  distance  of  3000  km.  The  shallow  depth  of  burial 


and  high  local  wave  velocity  (5.6  km/sec)  cause  P  and  pP  to 
arrive  almost  simultaneously  with  a  lag  of  =  0.17  seconds. 
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Fig.  4.41--First  mantle  arrival  for  Rio  Blanco.  Epicentral 
distance  3000  km,  source  S-l,  earth  structure  model  CIT  109P 
0.75  Q.  A.  =  62,  T1  =  0.6. 
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Fig.  4 . 42- -Reduced  velocity  potential  transform  amplitude  versus 
frequency  for  granite  source  functions  G-l  and  G-2,  scaled  to 
61  kt.  Also  shown  are  the  1  Hz  frequency  points  for  200,  600. 
and  1000  kt.  ’ 
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(a)  Source  G-l.  A^  -  42,  -  0.4,  A2  =  48,  T2  -  0.6. 
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(b)  Source  G-2.  A^  =  118,  =  0.4,  A2  ”  124,  T2  0.6. 


Fig.  4.43--First  mantle  arrival  synthetic  seismograms  for  Pile- 
driver  using  two  source  functions.  Yield  61  kt,  DOB  463  m, 
epicentral  distance  3000  km,  earth  structure  model  CIT  109P 
0.75  Q. 
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The  body  wave  magnitudes  for  this  event  are  summarized  in 
Table  4.3. 

The  magnitudes  measured  on  the  first  cycle,  m^^  ,  do 

not  relate  to  observations  and  are  given  only  for  comparison 

with  previous  sections.  The  synthetic  seismograms  are  quite 

similar  in  shape  to  the  first  two  cycles  of  the  actual  re- 
r 191  r 21 

cords1  1  and  m^  is  almost  certainly  directly  comparable 

to  the  LRSM  m,  . 

b 

As  was  the  case  for  the  sandstone  event  Rulison,  the 
maximum  amplitude  on  our  computed  records  (i,e.,  m£)  shows 
good  agreement  with  the  observed.  However,  once  again  the 
period  of  the  measured  cycle  is  somewhat  smaller  than  the 
periods  on  the  actual  records. 

Comparing  the  two  theoretical  records,  note  that  sub¬ 
stantial  differences  in  source  spectral  shape  have  little  in¬ 
fluence  on  the  shape  of  the  teleseismic  pulse.  The  combined 
filters  of  the  earth  structure  and  seismograph  have  too  nar¬ 
row  a  pass  band  to  see  much  source  spectral  variation.  The 
absolute  amplitude  of  the  source  spectrum  in  the  0.3  to  1  1 1  z 
range  is,  of  course,  directly  reflected  in  the  amplitude  of 
the  teleseismic  pulse. 

Magnitude  Versus  Yield  for  Granite 

In  Fig.  4.44,  synthetic  seismograms  for  explosions  of 
200,  600  and  1000  kt  in  granite  are  shown  along  with  a  plot 
of  magnitude  versus  yield.  The  changing  shape  of  the  source 
function  G-2  for  each  yield  is  indicated  on  Fig.  4.42.  All 
three  events  are  at  the  same  depth,  so  the  complicating  effect 
of  the  usual  functional  relationship  between  yield  and  DOB  is 
removed.  Despite  the  source  spectrum  differences,  the  slope 
of  the  resulting  magnitude  versus  log  yield  plot  (Fig.  4.44d) 
is  very  close  to  unity,  certainly  within  the  accuracy  of  our 
measurements . 
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(a)  200  kt.  A2  =  186,  T2  =  0.75. 


fl*  (ICCl 


(b)  600  kt.  =  564,  T2  =  0.75. 

Fig.  4. 44- -First  mantle  arrival  synthetic  seismograms  for  three 
explosions  in  granite.  Source  C-2,  DOB  1200  m,  epicentral 
distance  3000  km,  earth  structure  model  CIT  109P  0.75  Q. 
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The  Piled  river  event  was  characterized  by  a  substantial 
amount  of  anomalous  shear  radiation.  The  effect  of  a  rather 
large  tectonic  release  component  on  short  period  body  waves 
from  Pi led river  is  presented  in  Pig.  4.45  for  both  granite 
source  functions.  Assuming  the  tectonic  prestress  to  be 
aligned  along  a  north-south  axis,  as  is  typical  of  \'TS,  the 
chosen  azimuth  corresponds  to  stations  in  the  northeastern 
U.S.  and  adjacent  areas  of  Canada. 

I ven  for  this  large  tectonic  release  component  (the 
converted  sP^  is  of  the  same  order  as  direct  P  from  the 
explosion),  the  effect  on  SPZ  records  is  minimal.  The  mea¬ 
sured  magnitudes,  summarized  in  Table  4.9, are  seen  to  be 
little  different  than  the  corresponding  values  (Table  4.&) 
for  the  calculation  with  no  tectonic  release. 

4-4,5  Summary  of  Body  have  Magnitudes - -Comparison  to 

Observed 

The  body  wave  magnitudes  from  our  synthetic  seismo¬ 
grams  are  compared  to  observed  values  in  graphical  form  in 
Fig.  4.46.  The  computed  values  arc  all  from  computations 
using  the  earth  structure  model  CIT  109P  0.75  Q.  The  source 
function  T-4  was  used  for  all  tuff  calculations.  The  values 
for  Knickerbocker,  Greeley  and  Benhara  are  from  Table  4.6. 

The  Rulison  value  is  from  Table  4.7.  For  Piledriver  and 
Bilby  a  range  of  computed  values  is  indicated.  For  Pile- 
driver  this  represents  the  range  of  values  from  the  two 
source  functions  G-l  and  G-2,  with  and  without  tectonic 
release  (Tables  4.3  and  4.9).  For  Bilby  the  m^“^  values  arc 
taken  from  Table  4.5  which  gives  Bilby  mb  values  at  various 
teleseism ic  distances. 
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(a)  Source  G-l.  =  45,  -  0.4,  A2 


58,  T0  =  0.7. 
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(b)  Source  G-2.  A^  -  117, 


0.45,  A2  =  142,  T2  =  0.6. 


Fig.  4.45--First  mantle  arrival  synthetic  seismograms  for  the 
Piledriver  event  including  tectonic  release.  r;picentral  dis 
tance  3000  km,  aximuth  300°,  DOB  =  463  meters,  prestress  1 
-536  bars,  shatter  zone  radius  Rq  =  400  meters,  earth  struc¬ 
ture  model  CIT  109P  0.75  Q. 
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TABLH  4.9 

°!^Y  V.'AVn  MAGNITUDFS  FOR  PI  LFDR I  \'I:R  RECOUPS 
OF  FIGURF:  4.4  5 


Seismogram 

m,* 

b 

b 

(2) 
m,  J 
b 

4.45a 

(G-l) 

5.65 

5.91 

5.87 

4.4  5b 

(G-2) 

5.26 

5.55 

5.42 

"3  4-> 

<U 


4-1  ^ 
3  O 
GG 
E 

O  X) 
U  0) 
CO 

X)  3 

c 

rt  to 

3 

x  s: 

CD 

>  O' 

>H 

OJ  LO 

io  c- 

.o  • 

o  c 

*4-1  O, 
O  C3 

o 
c  —i 
o 

to  t— 

•r-(  h-H 

>1  u 

ca 

G  i — i 
E  <u 
O  "3 

u  o 

I  E 

i 

o  o 

*X  f-. 

•  3 

3-  4-i 

u 

•  3 

C/5  5- 


U-  CO 


153 


R- 2034 


V.  SUMMARY  AM!)  CONfLIISIOUS 

Using  a  deterministic  computer  model,  the  dependence 
of  teleseismic  ground  motion  on  the  near  source  material 
properties,  yield,  tectonic  stress  release,  depth  of  burial 
and  source - rece iver  earth  structure  has  been  investigated. 

In  Section  II  the  application  of  a  Lagrangian,  finite 
difference  shock  wave  code  (SKIPPER)  to  compute  the  equi¬ 
valent  elastic  source  representation  of  an  underground  ex¬ 
plosion  is  described.  Calculations  in  tuff,  sandstone  and 
granite  media  are  presented  demonstrating  the  effect  of 
changing  near-source  material  properties  on  the  explosion 
source  spectrum. 

In  Section  III  the  SKIPPER  code  output  in  the  form  of 
the  reduced  displacement  potential  O(t))  is  used  to  drive 
the  teleseismic  surface  wave  codes  to  compute  synthetic 
Rayleigh  wave  seismograms.  Determining  surface  wave  magni¬ 
tude,  Mg ,  directly  from  these  seismograms,  the  dependence 
of  this  important  parameter  on  the  source  spectrum,  depth  of 
burial  and  the  earth  structure  model  is  studied. 

Section  IV  is  devoted  to  a  study  of  the  dependence  of 
teleseismic  body  waves  on  the  important  controlling  para¬ 
meters.  The  explosion  source,  which  is  primarily  given  by 
the  SKIPPER  generated  reduced  displacement  potential,  may 
also  have  a  tectonic  stress  release  component.  Once  again 
synthetic  seismograms  are  generating  allowing  a  determina¬ 
tion  of  body  wave  magnitude,  m^ ,  in  a  standard  manner.  An 
investigation  of  the  dependence  of  m^  on  the  source  spectrum 
(including  tectonic  release),  yield,  depth  of  burial,  near 
source  wave  velocities  and  the  earth  structure  model  is 
presented . 

Among  the  important  conclusions  resulting  from  this 
research  are  the  following: 
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Civon  an  invariant  earth  structure  and  fixing  the 
depth  of  burial,  Rayleigh  wave  amplitude  is  linearly 
related  to  the  steady  state  value  of  the  reduced 
displacement  potential,  d(“). 

Rayleigh  wave  amplitude  is  nearly  insensitive  to 
changes  in  depth  of  burial  when  and  the  earth 

structure  are  invariant  for  DOB  in  the  range  of  0  to 
0400  ft. 

The  first  two  conclusions  imply  that  Rayleigh  wave 
amplitude  scales  as  the  yield,  V/,  of  the  explosion 
and  '1  scales  as  log  W  when  both  the  near  source 

•s  i  o 

material  properties  and  eartii  structure  are  in¬ 
var  i  an  t . 

Rayleigh  wave  magnitude  (M  )  is  quite  sensitive  to  the 
seismogram  phase  at  which  peak  to  peak  amplitude  is 
measured.  This  is  especially  true  if  the  earth  struc¬ 
ture  gives  a  large  Airy  phase  both  within  the  band 
pass  of  the  seismometer  and  near  20  seconds. 

For  the  first  arrival  P  wave  train,  the  spectral 
shape  and  position  of  the  spectral  peak  is  controlled 
almost  entirely  by  the  depth  of  burial.  The  wave 
train  is  assumed  to  be  composed  of  both  direct  and 
reflected  phases  and  the  shaping  of  the  spectrum 
is  caused  by  interference  between  the  two. 

Changing  the  explosion  yield  while  keeping  the  depth 
of  burial  fixed  does  not  appreciably  alter  the  shape 
of  the  P  wave  train.  In  this  situation  P  wave  ampli¬ 
tudes,  and  hence  m,  ,  scale  as  log  V.'. 

n  ’  l  o 

Tectonic  stress  release  at  the  source  should  have  little 
effect  on  m^.  The  major  tectonic  phase,  the  free  sur¬ 
face  converted  sP^,  phase,  arrives  too  late  to  influence 
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magnitude  measurements,  while  the  earlier  arriving 
PT  and  pPT  phases  are  generally  too  small  to 
significantly  alter  mb .  The  effect  of  tectonic 
stress  release  on  Mg  remains  an  area  for  future 
study . 

One  of  the  main  advantages  of  the  computational  model 
described  in  this  report  may  well  be  its  ability  to 
rapidly  invert  field  observations  to  deduce  the  earth 
structure  model  appropriate  for  a  specific  source- 
receiver  configuration.  Having  made  the  inversion 
to  obtain  reliable  earth  structure  models,  it  should 
then  be  possible  to  derive  the  explosion  yield  from 
teleseismic  data. 
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